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We consider the possibility of extracting the bottom quark mass from LEP data. 
The inclusive decay rate for Z — ► bb + ■ ■ ■ is obtained at order a s by summing up 
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the one-loop two-parton decay rate to the tree-level three-parton rate. We calcu- 
late the decay width of the Z-boson into two and three jets containing the 6-quark 
including complete quark mass effects. In particular, we give analytic results for 
a slight modification of the JADE clustering algorithm. We also study the angular 
distribution with respect to the angle formed between the gluon and the quark jets, 
which has a strong dependence on the quark mass. The impact of higher order QCD 
corrections on these observables is briefly discussed. Finally, we present numeri- 
cal results for some popular jet-clustering algorithms and show that, indeed, these 
three-jet observables are very sensitive to the 6-quark mass and well suited for its 
determination at LEP. 



CERN-TH.74 19/94 
February 1, 2008 



*) On leave of absence from the Joint Institute for Nuclear Research, Dubna, Russia. 
**) On leave of absence from Departament de Fisica Teorica, Universitat de Valencia, and IFIC, 
Valencia, Spain. 



1 Introduction 



In the Standard Model of electroweak interactions all fermion masses are free parameters and 
their origin, although linked to the spontaneous symmetry breaking mechanism, remains secret. 
Masses of charged leptons are well measured experimentally and neutrino masses, if they exist, 
are also bounded. In the case of quarks the situation is more complicated because free quarks 
are not observed in nature. Therefore, one can only get some indirect information on the values 
of the quark masses. For light quarks (m q < 1 GeV, the scale at which QCD interactions be- 
come strong), that is, for u-,d- and s-quarks, one can define the quark masses as the parameters 
of the Lagrangian that break explicitly the chiral symmetry of the massless QCD Lagrangian. 
Then, these masses can be extracted from a careful analysis of meson spectra and meson decay 
constants. For heavy quarks (c- and 6-quarks) one can obtain the quark masses from the known 
spectra of the hadronic bound states by using, e.g., QCD sum rules or lattice calculations. How- 
ever, since the strong gauge coupling constant is still large at the scale of heavy quark masses, 
these calculations are plagued by uncertainties and nonperturbative effects. 

It would be very interesting to have some experimental information on the quark masses 
obtained at much larger scales where a perturbative quark mass definition can be used and, 
presumably, non-perturbative effects are negligible. The measurements at LEP will combine 
this requirement with very high experimental statistics. 

The effects of quark masses can be neglected for many observables in LEP studies, as usu- 
ally quark masses appear in the ratio m 2 / m|. For the bottom quark, the heaviest quark produced 
at LEP, and taking a 6-quark mass of about 5 GeV this ratio is 0.003, even if the coefficient in 
front is 10 we get a correction of about 3%. Effects of this order are measurable at LEP, how- 
ever, as we will see later, in many cases the actual mass that should be used in the calculations 
is the running mass of the 6-quark computed at the mz scale: rhb(mz) ~ 3 GeV rendering the 
effect below the LEP precision for most of the observables. 

While this argument is correct for total cross sections for production of 6-quarks it is not 
completely true for quantities that depend on other variables. In particular it is not true for jet 
cross sections which depend on a new variable, y c (the jet-resolution parameter that defines the 
jet multiplicity) and which introduces a new scale in the analysis, E c = m z ^/y~ c - Then, for 
small values of y c there could be contributions coming like m\jE 2 = {m^jmz) 2 /y c which 
could enhance the mass effect considerably. In addition mass effects could also be enhanced by 
logarithms of the mass. For instance, the ratio of the phase space for two massive quarks and 
a gluon to the phase space for three massless particles is 1 + 8(m q /mz) 2 \og(m q /mz). This 
represents a 7% effect for m q = 5 GeV and a 3% effect for m q = 3 GeV. 

The high precision achieved at LEP makes these effects relevant. In fact, they have to be 
taken into account in the test of the flavour independence of a s (mz) 10-0] ■ In particular it has 
been shown [^] that the biggest systematic error in the measurement of a b s (mz) (ce s obtained 
from 66-production at LEP from the ratio of three to two jets) comes from the uncertainties in 
the estimate of the quark mass effects. This in turn means that mass effects have already been 
seen. Now one can reverse the question and ask about the possibility of measuring the mass of 
the bottom quark, mb, at LEP by assuming the flavour universality of the strong interactions. 

Such a measurement will also allow to check the running of m^ju) from fi = nib to /i = mz 
as has been done before for a s (/x). In addition fhb(mz) is the crucial input parameter in the 
analysis of the unification of Yukawa couplings predicted by many grand unified theories and 
which has attracted much attention in the last years [|7p . 

The importance of quark mass effects in Z-boson decays has already been discussed in the 
literature \^\. The complete order a s results for the inclusive decay rate of Z — > bb-\-bbg+- ■ ■ can 



i 



be found[] in [@]. The leading quark mass effects for the inclusive Z-width are known to order 
for the vector part [ |TT| ] and to order a 2 s for the axial- vector part [|12|]. Quark mass effects for 
three-jet final states in the process e + e~ — > qqg were considered first in [ |T3| ] for the photonic 
channel and extended later to the Z channel in [|14p and JT3|] . Recently [ |16| ] calculations of 
the three-jet event rates, including mass effects, were done for the most popular jet clustering 
algorithms using the Monte Carlo approach. 

In this paper we will discuss the possibility of measuring the 6-quark mass at LEP, in par- 
ticular, we study bottom quark mass effects in Z decays into two and three jets. In section 2 
we calculate the inclusive decay rate Z — > bb + bbg + • • • at order a s by summing one-loop 
virtual corrections to Z — > bb and the real gluon bremsstrahlung contribution. Dimensional 
continuation is used to regularize both infrared (IR) and ultraviolet (UV) divergences. Phase 
space integrations are also done in D dimensions. This calculation allows us to understand 
the details of the cancellation of IR divergences and how some, potentially large, logarithms 
of the quark mass are absorbed in the running quark mass fhb(mz)- In section 3 we calculate 
analytically the two and three-jet event rates in terms of the jet-resolution parameter y c and 
the mass of the quark for a slight modification of the well-known JADE algorithm Jl7| ] suit- 
able for analytic calculations with massive quarks. We also present numerical results for this 
scheme and for some of the most popular jet-clustering algorithms (DURHAM (K T ), JADE and 
E), estimate higher order contributions and compare with experimental results obtained by the 
DELPHI Collaboration [@] for 1990-1991 data. If the gluon jet can be identified with good effi- 
ciency a very interesting observable, which strongly depends on the quark mass, is the angular 
distribution with respect to the angle formed between the quark and the gluon jets. This distri- 
bution is calculated for massless quarks in section 4: analytically for JADE-type algorithms and 
numerically for the DURHAM algorithm. We also compute numerically the ratio of massive to 
massless angular distributions for the four jet-clustering algorithms. In section 5 we summarize 
the results obtained in the paper and comment on the possibility of using them to measure the 
6-quark mass in LEP experiments. Finally in the four appendices we collect all the functions 
and formulae needed in the body of the paper. 



2 The inclusive decay rate Z —> bb 

The main purpose of this paper is to investigate 6-quark mass effects in Z decays into two and 
three jets. Since at order a s the inclusive decay rate Z — > bb + ■ ■ ■ is given by the sum of the 
two- and three-jet decay widths we will start by studying this quantity. 

To calculate the total decay rate to order a s one has to sum up the virtual one-loop gluonic 
corrections to the Z — > bb with the real gluon bremsstrahlung. Both contributions are sepa- 
rately infrared divergent for massless gluons, therefore, some regularization method for the IR 
divergences is needed. The sum is, however, IR finite. 

Since there are many subtleties in this calculation, we sketch it in this section. Both pro- 
cesses, Z — > bb at one loop and Z — > bbg, are calculated in arbitrary dimension D = 4 — 2e and 
dimensional regularization is used to regularize the IR divergences [|18|]. At order a s and for 
massive quarks all IR divergences appear as simple poles 1/e. We show how the the divergences 
cancel in the sum and obtain the total inclusive rate. 

The first step is to compute the decay width Z — > bb at tree-level in dimension D. Since 
there are no IR divergences in this case it is not necessary to do the calculations in arbitrary 

1 The order a s corrections to the vector part, including the complete mass dependences, were already known 
from QED calculations jTo[|. 



space-time dimensions. However, there are IR divergences at the one-loop level and e factors 
could lead to finite contributions when multiplied by the divergent terms. 
The amplitude for the decay Z — > bb in D dimensions is 

T b (0) = f^u^^gy + g A l^{q) , (2.1) 

where the factor /i e has been included to make the gauge weak coupling g dimensionless in D 
dimensions; u x and v 2 are short-hand notations for the quark (antiquark) spinors, u\ = u(pi) 
and v 2 = v(p 2 ), e At (?) stands for the polarization vector of the Z-boson and gy (g A ) are the 
vector (axial-vector) neutral current couplings of the quarks in the Standard Model. At tree 
level and for the 6-quark we have 

4 2 

9v = -1 + g-V > 9a = 1 ■ (2.2) 

Here we denote by cw and sw the cosine and the sine of the weak mixing angle. 

Taking the square of the amplitude, averaging over initial state polarizations, summing over 
final state polarizations, and adding the phase space factor for the two-body decay given in 
appendix A QT~8| ] we obtain the following decay width in D dimensions, 



r[ 0) = OA/? 1 - 26 , (2.3) 



with 



C-m 92 r a-e) (miy 

and 

A b = - (3 2 - 2e)g 2 v + (3 2 (l - t)g\ . (2.5) 
In these expressions /3 is the relative velocity of the produced quarks 

2 

TYl 

f3 = VT^^n , n = -T- (2.6) 

m z z 

At the one-loop level (see diagrams in fig. [l]b), and after renormalization of the UV diver- 
gencesQ, the amplitude can be conveniently parameterized in terms of three form factors, f v , 
f A and fx, 

T b = (9v ((1 + \c 9 fv)l, + 4^ /t ^t) + + \c 9 fA)l,l,) v 2 e^q) , 

(2?) 

where C g is defined as follows, 

° 9 = v{£^) w^7)- (2 - 8) 

Here and below we will conventionally use a s = a s (mz) to denote the value of the running 
strong coupling at the m^-scale. 

2 Note that conserved currents or partially conserved currents as the vector and axial currents do not get renor- 
malized. Therefore, all UV divergences cancel when one sums properly self-energy and vertex diagrams. The 
remaining poles in e correspond to IR divergences. One can see this by separating carefully the poles correspond- 
ing to UV divergences from the poles corresponding to IR divergences. 
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The form factors, f v , /a and fr, are related by 

fv = fA + f T - (2.9) 

The two functions, fy and f A , contain an IR divergence, while f T is finite. Separating the 
divergent parts, we can rewrite the real parts of the form factors as follows (at order a s the 
imaginary parts will not contribute) 

Re{M = -~fe + fvf, (2-10) 

Re {/a} = ~fe + fAf, (2.11) 
Re{/ r } = f Tf , (2.12) 

where all functions f e , fvf, fAf and fxf are given in appendix B. Note that, as expected, the IR 
divergent part of the amplitude is proportional to the tree-level amplitude eq. fl2.ip. As the IR 
divergence manifests itself as a single pole in e, clearly, we only need to keep everywhere terms 
linear in e. 

From the amplitude ( |2~77| ) we obtain the one-loop corrected width in D dimensions 



with 



r D _ r (0) r (l) 
l b — l b + 1 b 



if > = -C g f t -Tf ] + C b C g (g v F v + g\F A ) , (2.13) 

where the finite functions F v and F A are given in appendix B in terms of the form factors and 
if ^ is given by eq. ( |23l ) 

The 0(a s ) result, eq. ( |2.13[ ), is divergent for e — > because the IR divergences associated 
with massless gluons running in the loops. To get a finite answer at this order we also need 
to include gluon bremsstrahlung from the quarks. This has to be computed by working in D 
dimensions. 

The amplitude for the process Z — > bbg (the two corresponding diagrams are given in fig. [j]c) 
can be written as 

T b) = » —g, „, y — 

+ U " v + ' JAl %X V + m>hV ) <(W(l) ■ (2-14) 

Here A a are the Gell-Mann SU (3) matrices, and e v a {k) is the gluon polarization vector. 

The square of the amplitude, in dimension D, gives a rather involved expression that can 
be conveniently simplified when one realizes that the most divergent part of it factorizes com- 
pletely, even in D dimensions, due to the factorization theorems for soft and collinear diver- 
gences. 

Adding the three-body phase space (see appendix A) we find that the decay width of Z — > 
bbg in D dimensions can be written as 

T bg = C b C g C F J d yi dy 2 6(h p )h-*A bg , (2.15) 

where C F = 4/3 is the SU(3) group factor, y x and y 2 are defined in terms of the energy fractions 
of the two outgoing quarks 

Vl = 2( Pl k)/m 2 z = 1 - 2E 2 /m z , y 2 = 2(p 2 k)/m 2 z = 1 - 2E l /m z (2.16) 



A 



and A bg comes from the square of the matrix element, 



h 



Ag = Ab-22 + 9v h V + 9A h A ■ 



(2.17) 



Here A b is the same combination of couplings and masses that appears in the tree-level decay 



width to two quarks, eq. (|2.5D, and the function h p is given by 



hp = 2/12/2(1 -yi- V2) - r b (yi + y-i)' 



(2.18) 



and it is exactly the same function that defines the phase space available for the three-body 
decay (see eq. ( [2.15| ) and appendix A). After phase space integration this term will contain an 
IR divergence which comes from the singularity at y\ = = 0. 

The functions hy and Ha describe the vector and the axial- vector parts of the remainder of 
the square of the amplitude which do not generate any IR divergence. In the limit e = they 
are given by: 



hv = o - + - ' (2 - 19) 

2 \yi yi) 

h A = (1 + 2r b ) h v + 2r b . (2.20) 
To perform the phase space integration it is convenient to change variables as follows 



with 



and 



yi 
2/2 

z - r b (l + z) 
z(l + z) 



g(z)w , 
g(z)zw , 



(z — c)(l — cz) 



(1 + c) 2 z(l + z) 



1 + p 



(2.21) 



(2.22) 



Then, both hy and Ha only depend on the variable z, and the function h p , which defines phase 
space and appears explicitly in eq. ( [2.17D , factorizes completely 



g{zfz(l + z)w 2 (l - w) . 



(2.23) 



The function g(z) has zeros at z\ = c and z 2 = 1/c. As phase space is defined by h p > we 
obtain that the phase space in terms of the new variables is given by 

c<z<l/c and < it; < 1 . (2.24) 

After this change of variables eq. Q2.15P can we rewritten as 

r 6g = C b C g C F [ 1/C dzg(z) 2 J 1 dwwh~ e A bg . (2.25) 

J c JO 

Now the w integration is very simple and leads to Beta functions. For the integration of the term 
of the amplitude proportional to A b (see eq. (|2.17|)) we get 



l/c f l 

dz I dwh 



p g(z) 2 z 2 w 3 



B(-2e,2-e) f dz^— Uzfz{\ + zjf * , (2.26) 
Jc q(z) z z z v > 



where the function £>(— 2e, 2 — e) has a single pole in e = 0. In this way, all the divergent 
behaviour has been factorized in the Beta function. Then, to perform the z integration we can 
expand the integrant for small e and keep only terms linear in e. The integrations can be easily 
performed and the results written in terms of logarithms and dilogarithmic functions. The rest 
of the integrals do not lead to any divergence and can be done, without problem, putting e equal 
to zero. 

After phase space integration, the decay width for Z — > bbg can be written in the following 
form 



1. 



bg 



Cgfe-^b + G b C g [g v G v + Q A G a 



(2.27) 



where the first term contains the IR divergent part and the IR finite functions Gy and Ga are 
given in appendix B. 

The IR divergent part of eq. ( |2.27D is identical, but with reversed sign, to the one obtained 
for , therefore in the sum they will cancel, as it should be: 



r (0) + r 



(i) 



r b9 = rf + C b C g {g 2 v T v + g 2 A T / 



(2.28) 



with 



Ty 



Fy + Gy , 

Fa + Ga- 



(2.29) 
(2.30) 



From the results of the appendix B we can easily obtain the limit of these functions for small 
quark masses, m b <C m z (fj, < 1) 



Ty 



1 + 12r 6 , 

1 - 6r fe (21ogr fe + 1] 



(2.31) 
(2.32) 



If we plug this result into eq. ( [2.28D we obtain the well-known result [ pT| ] 

,2 
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Z c^64tt 



gl (1 + ^(1 + 12r fe )j +g 2 A (l - Qr b + ^ (1 - Qr b (2 \ogr b + 1)) 



(2.33) 

It is interesting to note the presence of the large logarithm, log(m|/m^), proportional to the 
quark mass in the axial part of the QCD corrected width, eq. ( [2.33D . The mass that appears in 
all above calculations should be interpreted as the perturbative pole mass of the quark. But in 
principle the expression ( |2.33D could also be written in terms of the so-called running quark 
mass at the m z scale by using 



ml = m 2 b {m z ) 



1 + 2- 



.a.o 



7T 



log 



ml 



mi 



+ 



(2.34) 



Then, we see that all large logarithms are absorbed in the running of the quark mass from the 
m b scale to the m z scale [ITT]] and we have 



m z 



9 



c^64vr 



!Jv ( L + ^(1 + 12f 6 )J +9 A (l- Qf b + ^(1 



22n 



(2.35) 



where f b = ml(m z )/m 2 z . 

This result means that the bulk of the QCD corrections depending on the mass could be 
accounted for by using tree-level expressions for the decay width but interpreting the quark 



mass as the running mass at the mz scale. On the other hand, since fhb(mz) ~ 3 GeV is 
much smaller than the pole mass, m& ~ 5 GeV, it is clear that the quark mass corrections are 
much smaller than expected from the naive use of the tree-level result with m b m 5 GeV, which 
would give mass corrections at the 1.8% level while in fact, once QCD corrections are taken 
into account, the mass corrections are only at the 0.7% level. 

The final results of this section are well known but we find they could illuminate the dis- 
cussion of mass effects in the two- and three-jet event rates and in the angular distribution with 
respect to the angle formed between the quark and gluon jets. Moreover the intermediate results 
of this section will be used in the rest of the paper. 



3 Two- and three-jet event rates 

According to our current understanding of the strong interactions, coloured partons, produced 
in hard processes, are hadronized and, at experiment, one only observes colourless particles. It 
is known empirically that, in high energy collision, final particles group in several clusters by 
forming energetic jets, which are related to the primordial partons. Thus, in order to compare 
theoretical predictions with experiments, it is necessary to define precisely what is a jet in both, 
parton level calculations and experimental measurements. 

As we have seen in the previous section, at order a s , the decay widths of Z into both two and 
three partons are IR divergent. The two-parton decay rate is divergent due to the massless gluons 
running in the loops. The Z -boson decay width into three-partons has an IR divergence because 
massless gluons could be radiated with zero energy. The sum, however, is IR finite. Then it is 
clear that at the parton-level one can define an IR finite two-jet decay rate, by summing the two- 
parton decay rate and the IR divergent part of the three-parton decay width, e.g. integrated over 
the part of the phase space which contains soft gluon emission [|19[]. The integral over the rest 
of the phase space will give the three-jet decay rate. Thus we need to introduce a "resolution 
parameter" in the theoretical calculations in order to define IR-safe observables. Obviously, the 
resolution parameter, which defines the two- and the three-jet parts of the three-parton phase 
space should be related to the one used in the process of building jets from real particles. 

In the last years the most popular definitions of jets are based on the so-called jet clustering 
algorithms. These algorithms can be applied at the parton level in the theoretical calculations 
and also to the bunch of real particles observed at experiment. It has been shown that, for some 
of the algorithms, the passage from partons to hadrons (hadronization) does not change much 



the behaviour of the observables [pOQ, thus allowing to compare theoretical predictions with 
experimental results. In what follows we will use the word particles for both partons and real 
particles. 

In the jet-clustering algorithms jets are defined as follows: starting from a bunch of particles 
with momenta pj one computes, for example, a quantity like 

EE' 

Uij = 2—^(1 - cos%) 
s 

for all pairs (i, j) of particles. Then one takes the minimum of all and if it satisfies that it is 
smaller than a given quantity y c (the resolution parameter, y-cut) the two particles which define 
this tjij are regarded as belonging to the same jet, therefore, they are recombined into a new 
pseudoparticle by defining the four-momentum of the pseudoparticle according to some rule, 
for example 

Pk = Pi + Pj ■ 



Algorithm 


Resolution 


Combination 


EM 


2 (m)/ s 


Pk=Pi+ Pj 


JADE 


2{E i E J )/s (1 - cos%) 


Pk=Pi+ Pj 


E 


{Pi+Pjf/s 


Pk=Pi+ Pj 


DURHAM 


2min(X 2 , Ej)/s (1 - cos# -) 


Pk=Pi+ Pj 



Table 1 : The jet-clustering algorithms 

After this first step one has a bunch of pseudoparticles and the algorithm can be applied again 
and again until all the pseudoparticles satisfy > y c . The number of pseudoparticles found in 
the end is the number of jets in the event. 

Of course, with such a jet definition the number of jets found in an event and its whole 
topology will depend on the value of y c . For a given event, larger values of y c will result in a 
smaller number of jets. In theoretical calculations one can define cross sections or decay widths 
into jets as a function of y c , which are computed at the parton level, by following exactly the 
same algorithm. This procedure leads automatically to IR finite quantities because one excludes 
the regions of phase space that cause trouble. The success of the jet-clustering algorithms is due, 
mainly, to the fact that the cross sections obtained after the hadronization process agree quite 
well with the cross-sections calculated at the parton level when the same clustering algorithm 
is used in both theoretical predictions and experimental analyses. 

There are different successful jet-clustering algorithms and we refer to refs. [[201, j2l"l] for a 
detailed discussion and comparison of these algorithms in the case of massless quarks. 

In the rest of the paper we will use the four jet-clustering algorithms listed in the table |T|, 
where ^fs is the total centre of mass energy. In addition to the well-known JADE, E and 
DURHAM algorithms we will use a slight modification of the JADE scheme particularly useful 
for analytical calculations with massive quarks. It is defined by the two following equations 




Pk=Pi+ Pj 

We will denote this algorithm as the EM scheme. For massless particles and at the lowest 
order E, JADE and EM give the same answers. However already at order o? s they give different 
answers since after the first recombination the pseudoparticles are not massless anymore and 
the resolution functions are different. 

For massive quarks the three algorithms, E, JADE and EM are already different at order a s . 
The DURHAM (K T ) algorithm, which has been recently considered in order to avoid exponen- 
tiation problems present in the JADE algorithm []2(],^2p, is of course completely different from 
the other algorithms we use, both in the massive and the massless cases. 

In figure ^ we plotted the phase-space for two values of y c (y c = 0.04 and y c = 0. 14) for all 
four schemes (the solid line defines the whole phase space for Z — > qqg with m q = 10 GeV). 

There is an ongoing discussion on which is the best algorithm for jet clustering in the case 
of massless quarks. The main criteria followed to choose them are based in two requirements: 

1 . Minimize higher order corrections. 



2. Keep the equivalence between parton and hadronized cross sections. 

To our knowledge no complete comparative study of the jet-clustering algorithms has been done 
for the case of massive quarks. The properties of the different algorithms with respect to the 
above criteria can be quite different in the case of massive quarks from those in the massless 
case. The first one because the leading terms containing double-logarithms of y-cut (log 2 (?/ c )) 
that appear in the massless calculation (at order a s ) and somehow determine the size of higher 
order corrections are substituted in the case of massive quarks by single-logarithms of y c times 
a logarithm of the quark mass. The second one because hadronization corrections for massive 
quarks could be different from the ones for massless quarks. 

Therefore, we will not stick to any particular algorithm but rather present results and com- 
pare them for all the four algorithms listed in the table [j]. 



3.1 The analytic calculation for the EM scheme 

Here we calculate analytically, at leading order, the three-jet decay rate of the Z -boson by using 
the EM clustering algorithm. 

At the parton level the two-jet region in the decay Z — > bbg is given, in terms of the variables 
y 1 and y 2 , by the following conditions: 

yi<y c or y 2 < y c or l-2r h -y x -y 2 <y c . (3.1) 

This region contains the IR singularity, y\ — y 2 = and the rate obtained by the integration of 
the amplitude over this part of the phase space should be added to the one-loop corrected decay 
width for Z — > bb. The sum of these two quantities is of course IR finite and it is the so-called 
two-jet decay width at order a s . The integration over the rest of the phase space defines the 
three-jet decay width at the leading order. It is obvious that the sum of the two-jet and three-jet 
decay widths is independent of the resolution parameter y c , IR finite and given by the quantity 
T b = F(Z — ► bb + bbg + ■ • •) calculated in section 2. Therefore we have 

T b = T b 2j (y c ) + Tl J (y c ) + ... . 

Clearly, at order a s , knowing T b and T b 3j (y c ) we can obtain V b 2 -(y c ) as well. 

The calculation of Tlj(y c ) at order a s is a tree-level calculation and does not have any 
IR problem since the soft gluon region has been excluded from phase space. Therefore the 
calculation can be done in four dimensions without trouble. 

We will start with equation ( [2.1 5| ) taking e = and with the phase space constrained by the 



cuts defined in eq. flOj). 

1% = (c b C g C F J d yi dy 2 e PS e c A b ^ _ Q , (3.2) 

where the 9 function 

Ops = 0(h p ) (3.3) 
gives the whole phase space, and the product of 9 functions 

Oc = 0{y 2 - y c )6{ yi - y c )6{l - 2r b -y x -y 2 - y c ) (3.4) 

introduces the appropriate cuts for the EM scheme. The square of the amplitude, A bg , is given 
in eq. ( |2.17| ). The phase space and the cuts are represented in the first plot of fig. 



Q 



Depending on the value of y c the limits of integration are different, there are three cases 
which correspond to three different topologies of the overlapping of the phase space and the 
area defined by the cuts: 



Vc < 2r 6 
2r b < y c <y c 
Vc< Vc , (3.5) 

where y c = tJtHI — y^) + 0(r b 2 ^yr b ) is given by a solution of the following equation 

4(1 - 2y c - 2r b ) 2 (y 2 c + Ar b ) = y 2 c {2 -y B - 8r b ) 2 . (3.6) 

Since the integrant is symmetric under the exchange y x <-> y 2 we can restrict the region of 
integration to the region y x > y 2 (multiplying the result by a factor 2). In addition it is useful to 
change variables as before, y 2 = zy\. We will not discuss the technical details of the calculation 
here; all of the integrals can be reduced to logarithmic and dilogarithmic functions and the final 
result can be written in the following form 

Y% = C b C g (g 2 v H { °\y Cl r b )+g 2 A Hf(y Cl r b )) , (3.7) 

where the superscript (°) in the functions Hy^(y c , r b ) reminds us that this is only the lowest 
order result. Analytical expressions for the functions Hy \y c , r b ) and H^\y c , r b ) are given in 



appendix C. Obviously, the general form ( P.7| ) is independent of what particular jet-clustering 
algorithm has been used. 

In the limit of zero masses, r b = 0, chirality is conserved and the two functions Hy^ (y c , r b ) 
and (y c , r b ) become identical 

H$\y c ,0) = H%\y c ,0)=A(°\y c ). 

In this case we obtain the known result for the JADE-type algorithms, which is expressed in 
terms of the function (y c ) also given in appendix C Q. 

To see more clearly the size of mass effects we are going to study the following ratio of jet 
fractions 



p w _ ^j(yc)/^ b _ f H { y\y C) r b ) H^(y c ,r b ) i , , {] 



where we have defined 



9 V g 2 A 

° v ,2 _i_ „2 ' ° A „2 i „2 



9v + 9 a g v + g% 

In eq. ( pT8| ) we have kept only the lowest order terms in a s and r b . The last factor is due to 
the normalization to total rates. This normalization is important from the experimental point 
of view but also from the theoretical point of view because in these quantities large weak cor- 



rections dependent on the top quark mass 0230 cancel. Note that, for massless quarks, the ratio 



Tfj(y c )/T d is independent on the neutral current couplings of the quarks and, therefore, it is 
the same for up- and down-quarks and given by the function A^ . This means that we could 
equally use the normalization to any other light quark or to the sum of all of them (including 
also the c-quark if its mass can be neglected). 



3 Note that with our normalization A^(y c ) = hA(y c ), with A(y c ) defined in ref. [L 



in 



3.2 Estimate of higher order contributions 



All previous results come from a tree-level calculation, however, as commented in the introduc- 
tion, we do not know what is the value of the mass we should use in the final results since the 
difference among the pole mass, the running mass at [i = m b or the running mass at [i = mz 
are next-order effects in a s . 

In the case of the inclusive decay rate we have shown that one could account (with very 
good precision) for higher order corrections by using the running mass at the mz scale in the 
lowest order calculations. Numerically the effect of running the quark mass from m b to mz is 
very important. 

One could also follow a similar approach in the case of jet rates and try to account for 
the next-order corrections by using the running quark mass at different scales. We will see 
below that the dependence of R bd on the quark mass is quite strong (for all clustering schemes); 
using the different masses (e.g. mj or m b (Mz)) could amount to almost a factor 2 in the mass 
effect. This suggests that higher order corrections could be important. Here, however, the 
situation is quite different, since in the decay rates to jets we have an additional scale given by 
y c , E c = m z ^/y~ c , e.g. for y c = 0.01 we have E c = 9 GeV and for y c = 0.05, E c = 20 GeV. 
Perhaps one can absorb large logarithms, \og{m b /mz) by using the running coupling and the 
running mass at the fi = m z scale, but there will remain logarithms of the resolution parameter, 
log(y c ). For not very small y c one can expect that the tree-level results obtained by using 
the running mass at the mz scale are a good approximation, however, as we already said, the 
situation cannot be settled completely until a next-to-leading calculation including mass effects 
is available. 

Another way to estimate higher order effects in R\f is to use the known results for the 



massless case Q20, 21 , 24] 



Including higher order corrections the general form of eq. ( |3.7| ) is still valid with the change 



H v(A)(yci r b) -> H v(A) (y c ,r b ). Now we can expand the functions H v(A) (y c ,r b ) in a s and 
factorize the leading dependence on the quark mass as follows 

H v{ A)(y c ,r b ) = A^(y c ) + ^A^\y c ) + r b (4%(y c ,r 6 ) + ^B$ A) (y c ,r b )) +■■■ . (3.9) 

In this equation we already took into account that for massless quarks vector and axial contri- 
butions are identical^ 

Then, we can rewrite the ratio R%*, at order a s , as follows 



Rf = 1 + r b 



B ( y\y c ,r b ) (,a s ( 2#> (y c , r b ) A « (y c ) 



+CA 



AM(y c ) \ vr \B$\y c ,r b ) AW(y c ) 

B ( A\y c ,r b ) ( a s ( B A l \y C) r b ) A^{y c 
^\Vc) \ + vr \Bf(y c ,r b ) A<f»(y e 

x (i + 6r 6 (cA(l + 2^1og(r b ))-cv2^JJ . (3.10) 

From the calculations in this paper we know By\y c ,r b ) and B^\y c ,r b ); the lowest order 
function for the massless case, A^(y c ), is also known analytically for JADE-type algorithms, 

4 This is not completely true at 0(a 2 s ) because the triangle anomaly: there are one-loop triangle diagrams 
contributing to Z — > bbg with the top and the bottom quarks running in the loop. Since m t ^ mj, the anomaly 
cancellation is not complete. These diagrams contribute to the axial part even for = and lead to a deviation 
from Ay" 1 (y c ) — (y c ) [|5||. This deviation is, however, small [ ^5| ] and we are not going to consider its effect 
here. 
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eq. ( p,13|) and refs. [^D[|2T]], and for the DURHAM algorithm [JZ2p. A parameterization of the 
function (y c ) can be found in [ ^0[ ] for the different algorithms^. As we already mentioned 
this function is different for different clustering algorithms. The only unknown functions in 



and B A ^ (y c , r b ), which must be obtained from a complete calculation 



at order a 2 s including mass effects (at least at leading order in r b ). 

Nevertheless, in order to estimate the impact of higher order corrections in our calculation 
we will assume that B$ A {y c , r h )/B^ A {y c , r b ) < A^{y c )/A^{y c ) and take A^(y c )/A^(y c ) 
from[] [ ^0| , |2T| ]. Of course this does not need to be the case but at least it gives an idea of the 
size of higher order corrections. We will illustrate the numerical effect of these corrections for 
R\f in the next subsection. As we will see, the estimated effect of next-order corrections is 
quite large, therefore in order to obtain the 6-quark mass from these ratios the calculation of the 



functions By\(y c , r b ) is mandatory 



3.3 Numerical results for for different clustering algorithms 

To complete this section we present the numerical results for R\f calculated with the different 
jet-clustering algorithms. For the JADE, E and Durham algorithms we obtained the three-jet 
rate by a numerical integration over the phase-space given by the cuts (see fig. |2|). For the EM 
scheme we used our analytical results which were also employed to cross check the numerical 
procedure. 

In fig. ^| we present the ratio R\f, obtained by using the tree-level expression, eq. (|3T8|), 
against y c for m b = 5 GeV and m b = 3 GeV. We also plot the results given by eq. ( |3.10| ) 
(with By^ A (y c , r b )/By^ A (y c , r b ) = 0) for m b = 5 GeV, which gives an estimate of higher order 
corrections. For y c < 0.01 we do not expect the perturbative calculation to be valid. 

As we see from the figure, the behaviour of R\? is quite different in the different schemes. 
The mass effect has a negative sign for all schemes except for the E-algorithm. For y c > 0.05 
the mass effects are at the 4% level for m b = 5 GeV and at the 2% level for m b = 3 GeV (when 
the tree level expression is used). Our estimate of higher order effects, with the inclusion of the 
next-order effects in a s for massless quarks, shifts the curve for m b = 5 GeV in the direction of 
the 3 GeV result and amounts to about of 20% to 40% of the difference between the tree-level 
calculations with the two different masses. For both E and EM schemes we used the higher 
order results for the E scheme. 

For the JADE algorithm we have also plotted in fig. |3] the experimental results for R\f 
obtained by the DELPHI group [§] on the basis of the data collected in 1990-1991. The experi- 
mental errors, due to the limited statistics analyzed, are rather large. However, one can already 
see the effect of the quark mass. If the 6-quark mass would be zero, one should obtain a ratio 
R\? constant and equal to 1. It is clearly seen from the figure that for y c < 0.08 the data are 
significantly below 1. For larger values of y c , the number of events decreases, the errors become 
too large and the data are consistent with 1 . When larger amount of data is analyzed and the 
experimental error is decreased, it will be very interesting to see if data will exhibit the different 
signs of the mass effect in Rff (positive for the E scheme and negative for the other schemes) 
as predicted by our parton level calculations (see fig. |3|). 

In spite of the fact that the effect of the quark mass in has been seen, it is too early, in 
our opinion, to extract now the value of the 6-quark mass from the data. As discussed above 
the higher order corrections to R\f are presumably rather large and should be included in the 

5 With our choice of the normalization (y c ) = B(y c )/A, where B{y c ) is defined in [pop. 
6 For the EM algorithm this function has not yet been computed. To make an estimate of higher order corrections 
we will use in this case the results for the E algorithm. 
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Algorithm 


1.(0) 1.(2) 

Ky Ky hvy 


1.(0) l(1) l.(2) 
ft A "M ft A 


EM 


-2.72 -14.64 -28.58 


-2.61 -13.54 -30.67 


JADE 


-2.01 -5.19 -13.25 


-1.90 -4.13 -15.42 


E 


4.68 19.04 25.97 


4.71 19.81 23.39 


DURHAM 


-1.69 -4.76 -12.70 


-1.65 -4.28 -15.48 



Table 2: Results of the tree parameter fits of the functions By A (y c ,r b )/A^(y c ) = 
Elo lo S n 2/c in the ran S e °- 01 < Vc < 0-2 

theoretical calculations. However, it is clear, that once the essential next-to-leading order cor- 
rections will be available and all LEP data will be included in the analysis, the ratios R\f will 
certainly allow for a reasonable determination of the 6-quark mass and for a check of its running 
from m& to mz ■ 

To simplify the use of our results we present simple fits to the ratios By\(y c , r b ) /A^°\y c ), 
which define R^ d at lowest order, for the different clustering algorithms. We use the following 
parameterization: 

B^ A (y c ,r b )/A^(y c ) = £ k^\og n y c , (3.11) 

n=0 

and the results of the fits for the range 0.01 < y c < 0.2 are presented in table 0. 

In fig. |] we plot the ratios By A (y c , r b )/A^(y c ) as a function of y c for the different algo- 
rithms (dashed lines for m b = 5 GeV, dotted lines for m& = 3 GeV and solid curves for the 
result of our fits). As we see from the figure the remnant mass dependence in these ratios (in the 
range of masses we are interested in and in the range of y c we have considered) is rather small 
and for actual fits we used the average of the ratios for the two different masses. We see from 
these figures that such a simple three-parameter fit works reasonably well for all the algorithms. 

Concluding this section we would like to make the following remark. In this paper we 
discuss the Z-boson decay. In LEP experiments one studies the process e + e~ — > (Zj*) — > bb 
and, apart from the resonant Z-exchange cross section, there are contributions from the pure 
7-exchange and from the 7 — Z-interference. The non-resonant 7-exchange contribution at the 
peak is less than 1% for muon production and in the case of 6-quark production there is an 
additional suppression factor Q\ = 1/9. In the vicinity of the Z-peak the interference is also 
suppressed because it is proportional to Qh(s — m 2 z ) (y/s is the e + e~ centre of mass energy). We 
will neglect these terms as they give negligible contributions compared with the uncertainties in 
higher order QCD corrections to the quantities we are considering. 

Obviously, QED initial-state radiation should be taken into account in the real analysis; the 
cross section for 6-pair production at the Z resonance can be written as 

= / 4 b (s')F(s>/s)ds' (3.12) 



where F(s'/s) is the well-known QED radiator for the total cross section [ j2"7[ l and, the Born 



cross section, neglecting pure 7 exchange contribution and the 7 — Z-interference, has the form 

l2-KT e T b s 

abb[S) m 2 z (s-m|) 2 + m|r| K } 
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with obvious notation. Note that T b in this expression can be an inclusive width as well as some 
more exclusive quantity, which takes into account some kinematical restrictions on the final 
state. 



4 Angular distribution 

If the gluon jet can be identified with enough efficiency, an interesting quantity which is very 
sensitive to the IR behaviour of the amplitudes is the angular distribution with respect to the 
angle formed between one of the quark jets and the gluon jet[|. If -&i ($ 2 ) is the angle between 
the quark (antiquark) and the gluon jets we define d = min(i?i,i? 2 )- We want to obtain the 
angular distribution with respect to •&. The starting point is eq. ( p.2| ) where we change variables 
from one of the y\ or the y2 variables to To do this we take into account that the amplitude is 
completely symmetric in y\ and y 2 , therefore we can restrict the integration only to the region 
y 2 > yi and add a factor 2. In that region $ = i?x. Therefore to obtain the distribution with 
respect to -d it is enough to obtain the distribution with respect to "di but constraining the phase 
space integration to y 2 > y±. 

For y 2 > y\ we can easily express y\ in terms of cos •& = cos $1 as follows 



y 2 (l - y 2 - cos $J (1 - y 2 ) 2 - 4r fe ) 

y\ = — J =- L • (4.1) 

1 + 2/2 + cos d J (1 - y 2 ) 2 - Ar b 



Adding the Jacobian of the transformation we find from eq. fl3.2[ ) (taking e = as this quantity 
is IR convergent) 



dVl r 2y 2 J(l-y 2 ) 2 ~Ar b ) 

-f = C b C g C F 2 / dy 2 PS O c %2-</i) sintf- ^ J ^ A hg , (4.2) 

' ^l + 2/ 2 + cos^ v /(l-y 2 ) 2 -4V 



where y\ is expressed in terms of cos?? and y 2 using eq. (|4.1[). 

In order to see how large mass effects are in this angular distribution we define the following 
ratio of angular distributions: 

bd _ 1 dT tj I 1 dT % 

^ " p^r/ (4 - 3) 

In the case of massless quarks the integration limits in eq. Q4.2D can be found analytically 



for the JADE-type schemes and the result of the integration over y 2 is expressed in terms of 
logarithms involving ■d and y c . We find 



1 dY 



a, 



Y d d-d 7T 



U(Vc) , (4.4) 



where the function f^(y c ) is given analytically in appendix D for the JADE-type schemes and 
represented in fig. ^] for the JADE-type and the Durham algorithms for different values of y c 
(y c = 0.02 (solid line), y c = 0.04 (dashed line) y c = 0.06 (dotted line ) and y c = 0.08 
(dash-dotted line)). We observe a very sharp peak, for both algorithms, in the region of 90°- 
100° depending on the value of y c , for y c = 0.04 the peak is at about 92° for the JADE-type 
algorithms and at about 99° for the Durham algorithm. We see that the absolute size of the 
peak is a factor two larger in the case of the JADE-type algorithms (for the same value of y c ) 

7 We thank J. Fuster for suggesting us the study of this observable. 
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than in the case of the DURHAM scheme. This is due to the difference of phase spaces for two 
schemes. 

For massive quarks, although the integrations can still be performed analytically in the EM 
scheme, some of the integration limits are solutions of polynomial equations of the third degree 
and the analytical result is not especially enlighting. Then, we have computed the ratio R b $ by 
doing the one-dimensional integration in ( |4.2| ) numerically. 

Numerical results for R hd are presented in fig. |] for the different algorithms for y c = 0.04 
and for both m b = 5 GeV (solid line) and for m b = 3 GeV (dashed line). In all cases we plot 
the ratios for the interval of angles for which the differential cross section is still sizable (see 
fig. §), i.e. t? w 45° - 120° for JADE-type schemes and w 50° - 130° for the DURHAM 
algorithm. For small angles and m b = 5 GeV the effect can be as large as 10% of the ratio. 
Note, however, that the angular distribution, fig. |5], drops down rapidly for such small angles. In 
addition, since the ratio changes very fast in this region the exact size of the effect will depend 
on the angular resolution achieved at experiment. 

As in the case of ratios of three-jet event rates, R bd , the variation of the ratio of angular 
distributions, R$ d , for m b = 5 GeV and m b = 3 GeV gives a measure of the size of higher order 
corrections. 

We observe in all the ratios the irregular behaviour in the region where the massless angular 
distribution peaks. This is due to the fact that in the massive case the position of the peak is 
slightly shifted with respect to the massless case. The mismatch between the two peaks appears 
as a discontinuity in the ratio when seen from large scales. 

It will be interesting to see if data really follow these patterns for Rff. A preliminary analysis 
performed by the DELPHI group [ |28| ] seems to indicate that, indeed, data do follow these 
angular distributions, at least qualitatively, and exhibit the variations present in the different 
algorithms. 

5 Discussion and conclusions 

In this paper we have presented a theoretical study of quark-mass effects in the decay of the 
Z-boson into bottom quarks. 

First, we have reproduced, with the complete mass dependences, the results for the inclusive 
decay rate of the Z — > bb + ■ ■ ■ to order a s by adding gluon bremsstrahlung from the 6-quarks 
to the one-loop corrected decay width of Z — ► bb. Although the sum of the two contributions is 
finite, each of them is separately IR divergent. We used dimensional continuation to regularize 
the IR divergences and gave a complete analytical result in arbitrary space-time dimensions for 
each of the two contributions. 

The main contribution of this paper is, however, the analysis of some three-jet observables 
which are more sensitive to the value of the quark masses. 

For a slight modification of the JADE algorithm (the EM algorithm) we have calculated 
analytically the three-jet decay width of the Z-boson into 6-quarks as a function of the jet 
resolution parameter, y c , and the 6-quark mass. The answer is rather involved, but can be 
expressed in terms of elementary functions. Apart from the fact that these analytical calculations 
are interesting by themselves, they can also be used to test Monte Carlo simulations. For the 
EM, JADE, E and DURHAM clustering algorithms we have obtained the three-jet decay width 
by a simple two-dimensional numerical integration. Numerical and analytical results have been 
compared in the case of the EM scheme. 



We discussed quark-mass effects by considering the quantity 

which has many advantages from both the theoretical and the experimental point of views. In 
particular, at lowest order, the function F(m b , y c ) is almost independent on the quark mass (for 
the small values of the mass in which we are interested in) and has absolute values ranging from 
10 to 35 (depending on y c and on the algorithm), where the larger values are obtained for y c of 
about 0.01. 

At the lowest order in a s we do not know what is the exact value of the quark mass that 
should be used in the above equation since the difference between the different definitions of 
the 6-quark mass, the pole mass, m b « 5 GeV, or the running mass at the m^-scale, m b {m z ) ~ 
3 GeV, is order a s . Therefore, we have presented all results for these two values of the mass 
and have interpreted the difference as an estimate of higher order corrections. Conversely one 
can keep the mass fixed and include in F(m b , y c ) higher order corrections already known for 
the massless case. According to these estimates the 0(a s ) corrections can be about 40% of the 
tree-level mass effect (depending on the clustering scheme), although we cannot exclude even 
larger corrections. 

By using the lowest order result we find that for moderate values of the resolution parameter, 
y c ~ 0.05, the mass effect in the ratio R^ d is about 4% if the pole mass value of the 6-quark, 
m b « 5 GeV, is used, and the effect decreases to 2% if m b = 3 GeV. 

We have compared our predictions for R^ d for the JADE algorithm, with the results obtained 
from the 1990-1991 data by the DELPHI group [0. Although the errors obtained in the analysis 
of this limited sample of data are rather large, especially for y c > 0.08, one clearly sees that 
for small values of y-cut (y c < 0.08) the experimental points are systematically below 1, thus 
clearly exhibiting the effect of the mass of the quark, as for massless quark R^ d = 1. The 
size of the effect agrees roughly with the predictions. One can expect the reduction of the 
experimental error by, at least, a factor two when the data collected in 1992 are included in the 
analysis. Then, mass effects will be more clearly seen and it will be very interesting to see if 
data follow the different qualitative behaviour of the ratio R hd as a function of y c as predicted 
by the parton model calculations (positive effect for the E scheme and negative mass effect for 
the other algorithms). However, in order to extract a meaningful value of the fe-quark mass from 
the data it will be necessary to include next-to-leading order corrections since the leading mass 
effect we have calculated does not distinguish among the different definitions of the quark mass 
(pole mass, running mass at the m b scale or running mass at the rriz scale). We believe that the 
future analysis of the whole LEP statistics and its comparison with the theoretical predictions 
for the three-jet ratios, which meet the future experimental precision, will allow for a good 
determination of the 6-quark mass at the highest energy scale and for a check of its running 
from m b to mz- 

The high precision achieved at LEP allows for a good separation of the gluonic and quark 
jets and a measurement of the angular distribution of the radiated gluon with respect to the 
quark momenta. This angular distribution has been calculated for massless quarks analytically 
for the JADE-type schemes and numerically for the DURHAM algorithm. We have studied the 
mass effects, for the different jet-clustering algorithms, in the quantity 

1 dT b I 1 dT d 

r>bd _ 1 ai 3j / 1 ai 3j 



T b dd / T d dd 



We have shown that, for a reasonable value of the resolution parameter, y c = 0.04, the mass 
effects in this ratio can be as large as 10% of the ratio for m b = 5 GeV (depending on the 



algorithm, the angle d and the angular resolution). The larger values are obtained for small 
angles where, however, the angular distribution falls down very rapidly. A fit to this ratio can 
be used to extract the value of the 6-quark mass. It will be interesting to see if data really 
follow the predictions for the angular distributions and if the mass effects in the ratio of angular 
distributions are well described by our results. 

Concluding, we have raised the question of the possibility of measuring the 6-quark mass 
at LEP by using three-jet observables. In our opinion, this is a big challenge for both experi- 
mentalists and theorists. Clearly, more work has to be done in order the precision of theoretical 
predictions meet the experimental accuracy, in particular order a 2 s calculations and studies of 
hadronization corrections including mass effects will be needed. However, this effort is worth 
since it will allow for an independent measurement of m& at much larger energies where, pre- 
sumably, non-perturbative effects are negligible. 
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APPENDICES 



A Phase space in D = 4 — e dimensions 

The phase space for n-particles in the final state in D-dimensions [ |I"8l ] (D = 4 — 2e) has the 
following general form 

d(PS„) = (2n) D I] , C-fo F ^ (?- e p) ( a -d 



" (2tt)^2^ 



i=l, n 



{27v)D n (27r pi 2E / fa 2 - "flew** [?-Eft] ■ (A.2) 



i=l, n 



Then doing several trivial integrations we have the following phase-space factor for the process 

Z^bb 

PS - 1 i r(l-e) ^ 2 m|\- £ 

P52 -4^2r(2-2e) l^rj ' (A3) 

where /3 = \/I — 4r b with r fe = m^/m 2 z , 

For the case of the decay into three particles, Z — > we have 

where the function /i p which gives a phase-space boundary in terms of variables y\ = 2(pik)/m 2 z 
and y 2 = 2(p 2 k)/m 2 z has the form 

hp = 1/11/2(1 - V\ - 2/2) - n(yx + y 2 ) 2 ■ (A.5) 



B Inclusive decay rate functions 



In this section we collect the functions needed in section 2. The relevant form factors are: 

< 1 + /3 2 ^ 



U 
hs 

fAf 



2(3 



log(c) 



+ 



1 -3 2 

fJog(r b ) + C F 
l + P 2 < 



(B.l) 
(B.2) 



-2- 2 -±f MO 



7T 1 

Li 2 (c) + — - - log 2 (c) + log(c) log (1 - c) 



(B.3) 



In the expression for f A f, the first term, proportional to log(r b ), comes because our election for 
the term proportional to the divergence. The vector form factor, f V f can be written in terms of 
the other two form factors, 

fvj Im ■ In ■ (B.4) 
In terms of these form factors the functions Fy and Fa that appear in eq. Q2.13| ) are 



F v 



P'fAf ■ 



(B.5) 
(B.6) 
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The functions that come from real bremsstrahlung can be written as follows, 

G v = p{^(S-P 2 )G P + G Vh ^ , 



G A = (3 ((3 2 G P + G Ah ) , 



(B.7) 
(B.8) 



where 

G P = G Ph + 2/ e (l + log(/3)). (B.9) 

The terms proportional to f e come again from our choice of the coefficient of the divergence, 

and Gph is the finite part coming from the integration of the term proportional to h p in the 
amplitude 

G Ph = ^j c 1/C dzg(z)^\og(g(z) 3 z(l+z)) . (B.10) 
The result of the integration gives 



Gph — Cp 
1+/3 2 (I 



-2 log 



4/3 3 



l-(3 2 



+ 2-^log(c) 



; l 4 log 2 (c) + y - Li 2 (c) - Li 2 (c 2 ) - 3 log(c) log (1 + c) 



(B.ll) 



The functions G Vh and G Ah come from the integration of the h v and h A terms respectively 



9 - 2/3 2 + /? 4 
2^ 



log(c) , (B.12) 



and 



1 



where 



Ga/ 1 = §(3-/3 2 )G^+(1-/? 2 )G' a ,, 
^ f 1/c J / N2 / ^ 2 3 - 2/3 2 - /3 4 , 

Ah = 4^" L 9 ^ = T t 3 ~ ^ + log(c) 
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C Three-jet event rate functions 



The functions Hy* and H A °\ which give the leading contribution to the three-jet decay rate in 
the EM algorithm, can be written in the following form 
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H A °\y c ,n) 



Cf E ' 
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Cf E ' 
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K v + (1- (3 2 )K\ 



(C.l) 



with 



= 0(y c - y c ) 

02 = -0(y c -2r b )9(y c -y c ) 

03 = 0(y c - Vc) 



1 Q 



and y c ^Jrl(l — y/r^}. Here K l s corresponds to the soft part and K v ^ to the vector (axial) 
hard part. These functions are given by 



Kl = 
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where we used the following notation, 
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In the limit of massless quarks, r b — > 0, from the functions ify ■* and if^ ; given above we 
obtain 

r fe - 0) = H$\y c , r b ^0)^ A® , (C.12) 
Here the function (y c ) is the known result [ |Z(1| , |2"1| ] for the JADE algorithm 
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1 - Vc 
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(C.13) 



The function A(y c ) given in refs. [^0|,pTj differs from our A^°\y c ) in a factor 2 because we 
chose a different normalization for it. 



D Angular distribution functions 

The angular distribution studied in section 4 is given, in the massless case, by the function 
fti(y c )- In the JADE-type algorithms it can be written as follows 



i=l,2 



where the thetctj functions have the form 
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and 



h{Vc) = 



(i + by 
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(b + xA 



f2(y c ) = 

- 6(1 + 2b) log 
In these equations we defined 
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Figure 1 : Feynman diagrams contributing to the decay rates Z — > bb, Z — ► bbg at order a. 



Figure 2: The phase space for Z — ► 66(7 in the plane y\ and y 2 with cuts (y c = 0.04 and 
y c = 0.14) for the different algorithms. The mass of the quark has been set to 10 GeV to 
enhance mass effects in the plot. 



Figure 3: The ratios R^ d (see eq. (|3.8D) for the four algorithms. Solid lines correspond to 
m& = 5 GeV, dashed lines correspond to = 3 GeV and dotted lines give our estimate of 
higher order corrections to the m& = 5 GeV curve. For the JADE algorithm we have also 
included the results of the analysis of the data collected during 1990-1991 by the DELPHI 
group [§. 



Figure 4: The functions By /A^ and /A^ for the four algorithms. Dashed lines for 
m b = 3 GeV, dotted lines for m b = 5 GeV and solid lines for our three-parameter fit, eq. Q3.1 1|). 



Figure 5: Normalized angular distributions (eq. Q4.4[)) with respect to the angle formed between 
the quark and the gluon jets for the massless case for JADE-type and DURHAM algorithms. 
Solid line for y c = 0.02, dashed line for y c = 0.04, dotted line for y c = 0.06 and dash-dotted 
line for y c = 0.08 



Figure 6: The ratios of angular distributions R$ d (see eq. ([43|)) for y c = 0.04 for the different 



algorithms. Solid line for m b = 5 GeV and dashed line for m b = 3 GeV. 
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1 Introduction 



In the Standard Model of electroweak interactions all fermion masses are free parameters and 
their origin, although linked to the spontaneous symmetry breaking mechanism, remains secret. 
Masses of charged leptons are well measured experimentally and neutrino masses, if they exist, 
are also bounded. In the case of quarks the situation is more complicated because free quarks are 
not observed in nature. Therefore, one can only get some indirect information on the values of 
the quark masses. For light quarks (m q < 1 GeV, the scale at which QCD interactions become 
strong), that is, for u-,d- and s-quarks, one can define the quark masses as the parameters of 
the Lagrangian that break explicitly the chiral symmetry of the massless QCD Lagrangian. 
Then, these masses can be extracted from a careful analysis of meson spectra and meson decay 
constants. For heavy quarks (c- and 6-quarks) one can obtain the quark masses from the known 
spectra of the hadronic bound states by using, e.g., QCD sum rules or lattice calculations. 
However, since the strong gauge coupling constant is still large at the scale of heavy quark 
masses, these calculations are plagued by uncertainties and nonperturbative effects. 

It would be very interesting to have some experimental information on the quark masses 
obtained at much larger scales where a perturbative quark mass definition can be used and, 
presumably, non-perturbative effects are negligible. The measurements at LEP will combine 
this requirement with very high experimental statistics. 

The effects of quark masses can be neglected for many observables in LEP studies, as usually 
quark masses appear in the ratio m 2 /m 2 z . For the bottom quark, the heaviest quark produced at 
LEP, and taking a 6-quark mass of about 5 GeV this ratio is 0.003, even if the coefficient in front 
is 10 we get a correction of about 3%. Effects of this order are measurable at LEP, however, as 
we will see later, in many cases the actual mass that should be used in the calculations is the 
running mass of the 6-quark computed at the m z scale: fh b (m z ) ~ 3 GeV rendering the effect 
below the LEP precision for most of the observables. 

While this argument is correct for total cross sections for production of 6-quarks it is not 
completely true for quantities that depend on other variables. In particular it is not true for 
jet cross sections which depend on a new variable, y c (the jet-resolution parameter that defines 
the jet multiplicity) and which introduces a new scale in the analysis, E c = mz^fth- Then, 
for small values of y c there could be contributions coming like va\jE 2 = (m b /m z ) 2 /y c which 
could enhance the mass effect considerably. In addition mass effects could also be enhanced by 
logarithms of the mass. For instance, the ratio of the phase space for two massive quarks and 
a gluon to the phase space for three massless particles is 1 + 8(m q /mz) 2 log(m q /mz). This 
represents a 7% effect for m q = 5 GeV and a 3% effect for m q = 3 GeV. 

The high precision achieved at LEP makes these effects relevant. In fact, they have to be 
taken into account in the test of the flavour independence of a s (m z ) [1-5]. In particular it has 
been shown [6] that the biggest systematic error in the measurement of a b s (mz) (a s obtained 
from ^-production at LEP from the ratio of three to two jets) comes from the uncertainties in 
the estimate of the quark mass effects. This in turn means that mass effects have already been 
seen. Now one can reverse the question and ask about the possibility of measuring the mass of 
the bottom quark, m b , at LEP by assuming the flavour universality of the strong interactions. 

Such a measurement will also allow to check the running of rhb(fi) from fi = m,b to fi = mz 
as has been done before for a s (/i). In addition fhb(mz) is the crucial input parameter in the 
analysis of the unification of Yukawa couplings predicted by many grand unified theories and 
which has attracted much attention in the last years [7]. 

The importance of quark mass effects in Z-boson decays has already been discussed in the 
literature [8]. The complete order a s results for the inclusive decay rate of Z — >■ bb + bbg + • • • 



1 



can be found 1 in [9]. The leading quark mass effects for the inclusive zT-width are known to 
order a 3 s for the vector part [11] and to order a 2 s for the axial- vector part [12]. Quark mass effects 
for three-jet final states in the process e + e~ — >■ qqg were considered first in [13] for the photonic 
channel and extended later to the Z channel in [14] and [15]. Recently [16] calculations of 
the three-jet event rates, including mass effects, were done for the most popular jet clustering 
algorithms using the Monte Carlo approach. 

In this paper we will discuss the possibility of measuring the 6-quark mass at LEP, in 
particular, we study bottom quark mass effects in Z decays into two and three jets. In section 
2 we calculate the inclusive decay rate Z — >■ bb + bbg + • • • at order a s by summing one-loop 
virtual corrections to Z — >■ bb and the real gluon bremsstrahlung contribution. Dimensional 
continuation is used to regularize both infrared (IR) and ultraviolet (UV) divergences. Phase 
space integrations are also done in D dimensions. This calculation allows us to understand 
the details of the cancellation of IR divergences and how some, potentially large, logarithms 
of the quark mass are absorbed in the running quark mass m&(m#). In section 3 we calculate 
analytically the two and three-jet event rates in terms of the jet-resolution parameter y c and the 
mass of the quark for a slight modification of the well-known JADE algorithm [17] suitable for 
analytic calculations with massive quarks. We also present numerical results for this scheme 
and for some of the most popular jet-clustering algorithms (DURHAM (K T ), JADE and E), 
estimate higher order contributions and compare with experimental results obtained by the 
DELPHI Collaboration [2] for 1990-1991 data. If the gluon jet can be identified with good 
efficiency a very interesting observable, which strongly depends on the quark mass, is the 
angular distribution with respect to the angle formed between the quark and the gluon jets. 
This distribution is calculated for massless quarks in section 4: analytically for JADE-type 
algorithms and numerically for the DURHAM algorithm. We also compute numerically the 
ratio of massive to massless angular distributions for the four jet-clustering algorithms. In 
section 5 we summarize the results obtained in the paper and comment on the possibility of 
using them to measure the 6-quark mass in LEP experiments. Finally in the four appendices we 
collect all the functions and formulae needed in the body of the paper. 

2 The inclusive decay rate Z — >• bb 

The main purpose of this paper is to investigate 6-quark mass effects in Z decays into two and 
three jets. Since at order a s the inclusive decay rate Z — >■ bb + • • • is given by the sum of the 
two- and three-jet decay widths we will start by studying this quantity. 

To calculate the total decay rate to order a s one has to sum up the virtual one-loop gluonic 
corrections to the Z — >■ bb with the real gluon bremsstrahlung. Both contributions are sepa- 
rately infrared divergent for massless gluons, therefore, some regularization method for the IR 
divergences is needed. The sum is, however, IR finite. 

Since there are many subtleties in this calculation, we sketch it in this section. Both 
processes, Z — >■ bb at one loop and Z — >■ bbg, are calculated in arbitrary dimension D = 4 — 2e 
and dimensional regularization is used to regularize the IR divergences [18]. At order a s and for 
massive quarks all IR divergences appear as simple poles 1/e. We show how the the divergences 
cancel in the sum and obtain the total inclusive rate. 

The first step is to compute the decay width Z — >■ bb at tree-level in dimension D. Since 
there are no IR divergences in this case it is not necessary to do the calculations in arbitrary 

1 The order a s corrections to the vector part, including the complete mass dependences, were already known 
fromQED calculations [10]. 
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space-time dimensions. However, there are IR divergences at the one-loop level and e factors 
could lead to finite contributions when multiplied by the divergent terms. 
The amplitude for the decay Z — >■ bb in D dimensions is 

T b {0) = ^-r—^ii^av + gAH5)v 2 ^(q) , (2.1) 

where the factor // has been included to make the gauge weak coupling g dimensionless in D 
dimensions; u t and v 2 are short-hand notations for the quark (antiquark) spinors, u t = ■u(pi) 
and v 2 = v(p 2 ), ^{q) stands for the polarization vector of the zT-boson and gv (gA) are the 
vector (axial- vector) neutral current couplings of the quarks in the Standard Model. At tree level 
and for the 6-quark we have 

4 

9v = -1 + -s 2 w , gA = l ■ (2.2) 

Here we denote by cw and sw the cosine and the sine of the weak mixing angle. 

Taking the square of the amplitude, averaging over initial state polarizations, summing over 
final state polarizations, and adding the phase space factor for the two-body decay given in 
appendix A [18] we obtain the following decay width in D dimensions, 

r[ 0) = , (2.3) 

with 



- 2 r(l - e) (rn 
c 2 w MirT(2 - 2e) \4*fi 

and 



C b = m z ;r , (2.4) 



A b = 1(3 - (3 2 - 2e)g 2 v + - e)g\ . (2.5) 



In these expressions (3 is the relative velocity of the produced quarks 



(3 = yi^4^ , r b = ^ . (2.6) 



At the one-loop level (see diagrams in fig. lb), and after renormalization of the UV 
divergences 2 , the amplitude can be conveniently parameterized in terms of three form factors, 
fv, f A and/ r , 



T b = A* e ^«i (gv ((1 + \c 9 fv)l» + ^CJt^-) + g A (l + \c g f A )^ v 2 e»(q) , 

(2.7) 



4.c w V V 2 ' "* 2 yJ 2m b 
where C g is defined as follows, 




r(i - e ) 



(2.8) 



Here and below we will conventionally use a s = a s (m z ) to denote the value of the running 
strong coupling at the ra^-scale. 



2 Note that conserved currents or partially conserved currents as the vector and axial currents do not get 
renormalized. Therefore, all UV divergences cancel when one sums properly self-energy and vertex diagrams. The 
remaining poles in e correspond to IR divergences. One can see this by separating carefully the poles corresponding 
to UV divergences from the poles corresponding to IR divergences. 
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The form factors, f v , Ja and f T , are related by 



fv = fA + fr. (2.9) 

The two functions, f v and f A , contain an IR divergence, while fx is finite. Separating the 
divergent parts, we can rewrite the real parts of the form factors as follows (at order a s the 
imaginary parts will not contribute) 

Re{jV} = --Je + fvf, (2.10) 

Re{M = --Je + fAf, (2.11) 
Re{/r} = f Tf , (2.12) 

where all functions f e , fvf, fAf an d frf are given in appendix B. Note that, as expected, the IR 
divergent part of the amplitude is proportional to the tree-level amplitude eq. (2.1). As the IR 
divergence manifests itself as a single pole in e, clearly, we only need to keep everywhere terms 
linear in e. 

From the amplitude (2.7) we obtain the one-loop corrected width in D dimensions 



with 



rD _ r>(°) I rW 

L b — L b +!(, , 

lf } = -C g f e -T^ + C b C g (g 2 v F v + g 2 A F A ) , (2.13) 

where the finite functions Fy and Fa are given in appendix B in terms of the form factors and 
is given by eq. (2.3) 

The 0(a s ) result, eq. (2.13), is divergent for e — >■ because the IR divergences associated 
with massless gluons running in the loops. To get a finite answer at this order we also need 
to include gluon bremsstrahlung from the quarks. This has to be computed by working in D 
dimensions. 

The amplitude for the process Z — >■ bbg (the two corresponding diagrams are given in fig. lc) 
can be written as 



.2e 



9 - ( ¥ + m b )in(g v + g A j 5 ) 



+ M9r + 9il %X ¥ + mi)l " ) <(*KM • (2-i4) 

Here A a are the Gell-Mann SU(Z) matrices, and e v a {k) is the gluon polarization vector. 

The square of the amplitude, in dimension D, gives a rather involved expression that can be 
conveniently simplified when one realizes that the most divergent part of it factorizes completely, 
even in D dimensions, due to the factorization theorems for soft and collinear divergences. 

Adding the three-body phase space (see appendix A) we find that the decay width of Z — >■ bbg 
in D dimensions can be written as 

T bg = G b C g C F J d yi dy 2 6(h p )h;<A bg , (2.15) 

where Cf = 4/3 is the SU{2>) group factor, y x and y 2 are defined in terms of the energy fractions 
of the two outgoing quarks 

Vl = 2( Pl k)/m 2 z = 1 - 2E 2 /m z , y 2 = 2(p 2 k)/m 2 z = 1 - 2E 1 /m z (2.16) 
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and A bg comes from the square of the matrix element, 



A bg = A h -^ + g 2 v h v +g 2 A h A . (2. 17) 

2/12/2 

Here A b is the same combination of couplings and masses that appears in the tree-level decay 
width to two quarks, eq. (2.5), and the function h p is given by 

K = 2/12/2(1 - yi - 2/2) - n(yi + 2/2) 2 , (2.18) 

and it is exactly the same function that defines the phase space available for the three-body 
decay (see eq. (2.15) and appendix A). After phase space integration this term will contain an 
IR divergence which comes from the singularity at y t = y 2 = 0. 

The functions h v and h A describe the vector and the axial- vector parts of the remainder of 
the square of the amplitude which do not generate any IR divergence. In the limit e = they 
are given by: 

kv = + , (2.19) 

2 \yi 2/2/ 

h A = (1 + 2r b )h v + 2r b . (2.20) 
To perform the phase space integration it is convenient to change variables as follows 

2/1 = g(z)w , 

2/2 = g{z)zw , 

with 

g(z) = £^(i±£)! = 1 (£ - - <*) (2.21) 

9{ ' z(l + z) (1 + c) 2 z(l + z) K ) 

and 

c = ^ . (2.22) 

Then, both h v and h A only depend on the variable z, and the function h p , which defines phase 
space and appears explicitly in eq. (2.17), factorizes completely 

h p = g{zfz{l + z)w 2 (l - w) . (2.23) 

The function g(z) has zeros at z\ = c and z 2 = 1/c. As phase space is defined by h p > we 
obtain that the phase space in terms of the new variables is given by 

c < z < 1/c and < w < 1 . (2.24) 

After this change of variables eq. (2.15) can we rewritten as 

T bg = C b C g C F dzg(z) 2 £ dwwh; e A bg . (2.25) 

Now the w integration is very simple and leads to Beta functions. For the integration of the term 
of the amplitude proportional to A b (see eq. (2.17)) we get 
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where the function B(—2e,2 — e) has a single pole in e = 0. In this way, all the divergent 
behaviour has been factorized in the Beta function. Then, to perform the z integration we can 
expand the integrant for small e and keep only terms linear in e. The integrations can be easily 
performed and the results written in terms of logarithms and dilogarithmic functions. The rest 
of the integrals do not lead to any divergence and can be done, without problem, putting e equal 
to zero. 

After phase space integration, the decay width for Z — >■ bbg can be written in the following 
form 

T bg = C g f e -V { b 0) + C b C g (g 2 v G v + g 2 A G A ) , (2.27) 

where the first term contains the IR divergent part and the IR finite functions Gv and Ga are 
given in appendix B. 

The IR divergent part of eq. (2.27) is identical, but with reversed sign, to the one obtained 
for therefore in the sum they will cancel, as it should be: 



r (o) , r (i) , r 



if } + C b C g (g 2 v T v + g 2 A T A ) , 



(2.28) 



with 



Fy + Gv 
Fa + Ga 



(2.29) 
(2.30) 



From the results of the appendix B we can easily obtain the limit of these functions for small 
quark masses, m b <C m z (r b <C 1) 



1 + 12r 6 , 

l-6r 6 (21og r 6 + l) 



(2.31) 
(2.32) 



If we plug this result into eq. (2.28) we obtain the well-known result [11] 

„2 



m z 



c 2 w Miv 



g 2 v (l + ^(1 + 12r 6 )) +g 2 A (l - 6r b + ^ (1 - 6r 6 (21og r b + 1))) . 

(2.33) 

It is interesting to note the presence of the large logarithm, log(m^/m|), proportional to the 
quark mass in the axial part of the QCD corrected width, eq. (2.33). The mass that appears in 
all above calculations should be interpreted as the perturbative pole mass of the quark. But in 
principle the expression (2.33) could also be written in terms of the so-called running quark 
mass at the m z scale by using 



m\ = m 2 b (m z ) 



a 



1 + 2- log -f +o 



7T 



mi 



(2.34) 



Then, we see that all large logarithms are absorbed in the running of the quark mass from the 
m b scale to the m z scale [11] and we have 



m z 



c 2 w Mir 



g 2 v (l + ^(1 + 12r 6 )) + g\ (l " 6r & + ^(1 - 22r 6 )) 



(2.35) 



where f b = ml(m z )/m 2 z . 

This result means that the bulk of the QCD corrections depending on the mass could be 
accounted for by using tree-level expressions for the decay width but interpreting the quark 
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mass as the running mass at the m z scale. On the other hand, since rh b {m z ) ~ 3 GeV is 
much smaller than the pole mass, m b 5 GeV, it is clear that the quark mass corrections are 
much smaller than expected from the naive use of the tree-level result with m b 5 GeV, which 
would give mass corrections at the 1.8% level while in fact, once QCD corrections are taken 
into account, the mass corrections are only at the 0.7% level. 

The final results of this section are well known but we find they could illuminate the 
discussion of mass effects in the two- and three-jet event rates and in the angular distribution 
with respect to the angle formed between the quark and gluon jets. Moreover the intermediate 
results of this section will be used in the rest of the paper. 

3 Two- and three-jet event rates 

According to our current understanding of the strong interactions, coloured partons, produced 
in hard processes, are hadronized and, at experiment, one only observes colourless particles. It 
is known empirically that, in high energy collision, final particles group in several clusters by 
forming energetic jets, which are related to the primordial partons. Thus, in order to compare 
theoretical predictions with experiments, it is necessary to define precisely what is a jet in both, 
parton level calculations and experimental measurements. 

As we have seen in the previous section, at order a s , the decay widths of Z into both two and 
three partons are IR divergent. The two-parton decay rate is divergent due to the massless gluons 
running in the loops. The zT-boson decay width into three-partons has an IR divergence because 
massless gluons could be radiated with zero energy. The sum, however, is IR finite. Then it 
is clear that at the parton-level one can define an IR finite two-jet decay rate, by summing the 
two-parton decay rate and the IR divergent part of the three-parton decay width, e.g. integrated 
over the part of the phase space which contains soft gluon emission [19]. The integral over the 
rest of the phase space will give the three-jet decay rate. Thus we need to introduce a "resolution 
parameter" in the theoretical calculations in order to define IR-safe observables. Obviously, the 
resolution parameter, which defines the two- and the three-jet parts of the three-parton phase 
space should be related to the one used in the process of building jets from real particles. 

In the last years the most popular definitions of jets are based on the so-called jet clustering 
algorithms. These algorithms can be applied at the parton level in the theoretical calculations 
and also to the bunch of real particles observed at experiment. It has been shown that, for some 
of the algorithms, the passage from partons to hadrons (hadronization) does not change much 
the behaviour of the observables [20], thus allowing to compare theoretical predictions with 
experimental results. In what follows we will use the word particles for both partons and real 
particles. 

In the jet-clustering algorithms jets are defined as follows: starting from a bunch of particles 
with momenta pi one computes, for example, a quantity like 

Vij = 2 (1 - cos 6ij) 

s 

for all pairs (z, j) of particles. Then one takes the minimum of all y^ and if it satisfies that it is 
smaller than a given quantity y c (the resolution parameter, y-cut) the two particles which define 
this yij are regarded as belonging to the same jet, therefore, they are recombined into a new 
pseudoparticle by defining the four-momentum of the pseudoparticle according to some rule, 
for example 

Pk = Pi + Pj ■ 
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Algorithm 


Resolution 


Combination 


EM 


2 {PiPj)/ s 


Pk = Pi + Pj 


JADE 


2(E i E j )/s(l-cos#i j ) 


Pk=Pi+ Pj 


E 


(Pi+Pj)V s 


Pk=Pi+ Pj 


DURHAM 


2mm(E?,E])/s(l- cos tfy) 


Pk=Pi+ Pj 



Table 1: The jet-clustering algorithms 

After this first step one has a bunch of pseudoparticles and the algorithm can be applied again 
and again until all the pseudoparticles satisfy yij > y c . The number of pseudoparticles found in 
the end is the number of jets in the event. 

Of course, with such a jet definition the number of jets found in an event and its whole 
topology will depend on the value of y c . For a given event, larger values of y c will result in a 
smaller number of jets. In theoretical calculations one can define cross sections or decay widths 
into jets as a function of y c , which are computed at the parton level, by following exactly the 
same algorithm. This procedure leads automatically to IR finite quantities because one excludes 
the regions of phase space that cause trouble. The success of the jet-clustering algorithms is due, 
mainly, to the fact that the cross sections obtained after the hadronization process agree quite 
well with the cross-sections calculated at the parton level when the same clustering algorithm is 
used in both theoretical predictions and experimental analyses. 

There are different successful jet-clustering algorithms and we refer to refs. [20,21] for a 
detailed discussion and comparison of these algorithms in the case of massless quarks. 

In the rest of the paper we will use the four jet-clustering algorithms listed in the table 1, 
where y/s is the total centre of mass energy. In addition to the well-known JADE, E and 
DURHAM algorithms we will use a slight modification of the JADE scheme particularly useful 
for analytical calculations with massive quarks. It is defined by the two following equations 

y v = 2^ 

S 

and 

Pk = Pi + Pj 

We will denote this algorithm as the EM scheme. For massless particles and at the lowest 
order E, JADE and EM give the same answers. However already at order a 2 s they give different 
answers since after the first recombination the pseudoparticles are not massless anymore and 
the resolution functions are different. 

For massive quarks the three algorithms, E, JADE and EM are already different at order 
a s . The DURHAM (Kt) algorithm, which has been recently considered in order to avoid 
exponentiation problems present in the JADE algorithm [22,20] , is of course completely different 
from the other algorithms we use, both in the massive and the massless cases. 

In figure 2 we plotted the phase-space for two values of y c (y c = 0.04 and y c = 0.14) for all 
four schemes (the solid line defines the whole phase space for Z — >■ qqg with m q = 10 GeV). 

There is an ongoing discussion on which is the best algorithm for jet clustering in the case 
of massless quarks. The main criteria followed to choose them are based in two requirements: 

1 . Minimize higher order corrections. 
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2. Keep the equivalence between parton and hadronized cross sections. 

To our knowledge no complete comparative study of the jet-clustering algorithms has been done 
for the case of massive quarks. The properties of the different algorithms with respect to the 
above criteria can be quite different in the case of massive quarks from those in the massless 
case. The first one because the leading terms containing double-logarithms of y-cut (log 2 (y c )) 
that appear in the massless calculation (at order a s ) and somehow determine the size of higher 
order corrections are substituted in the case of massive quarks by single-logarithms of y c times 
a logarithm of the quark mass. The second one because hadronization corrections for massive 
quarks could be different from the ones for massless quarks. 

Therefore, we will not stick to any particular algorithm but rather present results and compare 
them for all the four algorithms listed in the table 1 . 

3.1 The analytic calculation for the EM scheme 

Here we calculate analytically, at leading order, the three-jet decay rate of the zT-boson by using 
the EM clustering algorithm. 

At the parton level the two-jet region in the decay Z — >■ bbg is given, in terms of the variables 
yi and y 2 , by the following conditions: 

yi < y c or y 2 < y c or 1 - 2r b - y x - y 2 < y c . (3.1) 

This region contains the IR singularity, y t = y 2 = and the rate obtained by the integration of 
the amplitude over this part of the phase space should be added to the one-loop corrected decay 
width for Z — >■ bb. The sum of these two quantities is of course IR finite and it is the so-called 
two-jet decay width at order a s . The integration over the rest of the phase space defines the 
three-jet decay width at the leading order. It is obvious that the sum of the two-jet and three-jet 
decay widths is independent of the resolution parameter y c , IR finite and given by the quantity 
Tb = T(Z — >■ bb + bbg + • • •) calculated in section 2. Therefore we have 

T b = T b 2j (y c ) + T b 3j (y c ) + ■■■ . 

Clearly, at order a s , knowing T b and Y\-{y c ) we can obtain Y b 2 -{y c ) as well. 

The calculation of Y\-{y c ) at order a s is a tree-level calculation and does not have any 
IR problem since the soft gluon region has been excluded from phase space. Therefore the 
calculation can be done in four dimensions without trouble. 

We will start with equation (2.15) taking e = and with the phase space constrained by the 
cuts defined in eq. (3.1). 

1% = (c b C g C F J d yi dy 2 e PS e c A h ^ _ Q , (3.2) 

where the 9 function 

Ops = 0(h p ) (3.3) 
gives the whole phase space, and the product of 6 functions 

Oc = 0(y 2 - y c )6( yi - y c )6(l - 2r b - y x - y 2 - y c ) (3.4) 

introduces the appropriate cuts for the EM scheme. The square of the amplitude, A bg , is given 
in eq. (2.17). The phase space and the cuts are represented in the first plot of fig. 2. 
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Depending on the value of y c the limits of integration are different, there are three cases 
which correspond to three different topologies of the overlapping of the phase space and the 
area defined by the cuts: 

y c < 2r 6 
2r 6 < y c < y c 
Vc< Vc , (3.5) 

where y c = y/rb(l — + 0(rl is given by a solution of the following equation 

4(1 - 2y c - 2r h )\y 2 c + 4r 6 ) = y 2 c (2 - y c - 8r b ) 2 . (3.6) 

Since the integrant is symmetric under the exchange y 1 <->■ y 2 we can restrict the region of 
integration to the region y t > y 2 (multiplying the result by a factor 2). In addition it is useful to 
change variables as before, y 2 = zy t . We will not discuss the technical details of the calculation 
here; all of the integrals can be reduced to logarithmic and dilogarithmic functions and the final 
result can be written in the following form 

1% = C b C g (g 2 v H$\y c , r b ) + g 2 A H A °\y c , r b j) , (3.7) 

where the superscript (°) in the functions H^ A -^(y c , r b ) reminds us that this is only the lowest 

order result. Analytical expressions for the functions H$\y c , r b ) and H A °\y c , r b ) are given in 
appendix C. Obviously, the general form (3.7) is independent of what particular jet-clustering 
algorithm has been used. 

In the limit of zero masses, r b = 0, chirality is conserved and the two functions H^\y c , r b ) 
and H A °\y c , r b ) become identical 

H$\y e ,0) = H A °\y c ,0) = A<°>(y c ) . 

In this case we obtain the known result for the JADE-type algorithms, which is expressed in 
terms of the function A^°\y c ) also given in appendix C 3 . 

To see more clearly the size of mass effects we are going to study the following ratio of jet 
fractions 

R ' = fJfcjT^ = [ CV ^M- + CA ^M~ j ( + 1 A + ( b)] ' ' ' 

where we have defined 

V ~ 2 , 2 ' A ~ 2 , 2 - 

9v +9a 9v+ 9 a 

In eq. (3.8) we have kept only the lowest order terms in a s and r b . The last factor is due to the 
normalization to total rates. This normalization is important from the experimental point of view 
but also from the theoretical point of view because in these quantities large weak corrections 
dependent on the top quark mass [23] cancel. Note that, for massless quarks, the ratio T^- (y c ) /T d 
is independent on the neutral current couplings of the quarks and, therefore, it is the same for 
up- and down-quarks and given by the function A^°\ This means that we could equally use the 
normalization to any other light quark or to the sum of all of them (including also the c-quark if 
its mass can be neglected). 



3 Note that with our normalization A^°\y c ) = \A{y c ), with A{y c ) defined in ref. [20]. 
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3.2 Estimate of higher order contributions 



All previous results come from a tree-level calculation, however, as commented in the introduc- 
tion, we do not know what is the value of the mass we should use in the final results since the 
difference among the pole mass, the running mass at /j, = m b or the running mass at /j, = m z 
are next-order effects in a s . 

In the case of the inclusive decay rate we have shown that one could account (with very 
good precision) for higher order corrections by using the running mass at the mz scale in the 
lowest order calculations. Numerically the effect of running the quark mass from m b to mz is 
very important. 

One could also follow a similar approach in the case of jet rates and try to account for the 
next-order corrections by using the running quark mass at different scales. We will see below that 
the dependence of R bd on the quark mass is quite strong (for all clustering schemes); using the 
different masses (e.g. m b or m b (M z )) could amount to almost a factor 2 in the mass effect. This 
suggests that higher order corrections could be important. Here, however, the situation is quite 
different, since in the decay rates to jets we have an additional scale given by y c , E c = m Z yJy~ c , 
e.g. for y c = 0.01 we have E c = 9 GeV and for y c = 0.05, E c = 20 GeV. Perhaps one can 
absorb large logarithms, log(m b /mz) by using the running coupling and the running mass at 
the /j, = mz scale, but there will remain logarithms of the resolution parameter, log(y c ). For 
not very small y c one can expect that the tree-level results obtained by using the running mass 
at the mz scale are a good approximation, however, as we already said, the situation cannot be 
settled completely until a next-to-leading calculation including mass effects is available. 

Another way to estimate higher order effects in R hd is to use the known results for the 
massless case [24,20,21] 

Including higher order corrections the general form of eq. (3.7) is still valid with the change 
H^ A j(y c ,r b ) ->■ H V ( A -)(y c ,r b ). Now we can expand the functions H V ( A )(y c , r b ) in a s and 
factorize the leading dependence on the quark mass as follows 

H V (A)(yc,r b ) = A<°>(y c ) + ^-A^(y c ) + r b {B^ {A) {y c ,r b ) + ^B$ A) (y c ,r b ^ +■■■ . (3.9) 

In this equation we already took into account that for massless quarks vector and axial contri- 
butions are identical 4 

Then, we can rewrite the ratio R b 3 d , at order a s , as follows 

B A °\y c ,r b ) ( a s ( B<}\y e ,r b ) AW(y c ) \\ 
+CA A(o)(y c ) ^ + , \B A % c ,r b ) A(o)(y c )))_ 

x (l + 6r b (c A {l + 2^1og(r 6 )) - c v 2^j . (3.10) 

From the calculations in this paper we know By\y c , r b ) and B A °\y c , r b ); the lowest order 
function for the massless case, A^°\y c ), is also known analytically for JADE-type algorithms, 
eq. (C.13) and refs. [20,21], and for the DURHAM algorithm [22]. A parameterization of the 

4 This is not completely true at 0{a 2 s ) because the triangle anomaly: there are one-loop triangle diagrams 
contributing to Z — > bbg with the top and the bottom quarks running in the loop. Since m t ^ mj, the anomaly 
cancellation is not complete. These diagrams contribute to the axial part even for mj, = and lead to a deviation 
from Ay\y c ) = A A \y c ) [25]. This deviation is, however, small [25] and we are not going to consider its effect 
here. 



11 



function A^\y c ) can be found in [20] for the different algorithms 5 . As we already mentioned 
this function is different for different clustering algorithms. The only unknown functions in 
eq. (3.10) are By\y c , r b ) and B A \y c ,r b ), which must be obtained from a complete calculation 
at order a 2 s including mass effects (at least at leading order in r b ). 

Nevertheless, in order to estimate the impact of higher order corrections in our calculation 
we will assume that B^ A (y c ,r b )/B^ A (y c ,r b ) < A^\y c )l A^\y c ) and take A^\y c )l A^\y c ) 
from 6 [20, 21]. Of course this does not need to be the case but at least it gives an idea of the 
size of higher order corrections. We will illustrate the numerical effect of these corrections for 
R bd in the next subsection. As we will see, the estimated effect of next-order corrections is 
quite large, therefore in order to obtain the 6-quark mass from these ratios the calculation of the 
functions By^ A (y c ,r b ) is mandatory [26]. 

3.3 Numerical results for R b / for different clustering algorithms 

To complete this section we present the numerical results for R 3 d calculated with the different 
jet-clustering algorithms. For the JADE, E and Durham algorithms we obtained the three-jet 
rate by a numerical integration over the phase-space given by the cuts (see fig. 2). For the EM 
scheme we used our analytical results which were also employed to cross check the numerical 
procedure. 

In fig. 3 we present the ratio R b 3 d , obtained by using the tree-level expression, eq. (3.8), 
against y c for m b = 5 GeV and m b = 3 GeV. We also plot the results given by eq. (3.10) 
(with B^ A (y c , r b )/ B$ A (y c , r b ) = 0) for m b = 5 GeV, which gives an estimate of higher order 
corrections. For y c < 0.01 we do not expect the perturbative calculation to be valid. 

As we see from the figure, the behaviour of R b 3 d is quite different in the different schemes. 
The mass effect has a negative sign for all schemes except for the E-algorithm. For y c > 0.05 
the mass effects are at the 4% level for m b = 5 GeV and at the 2% level for m b = 3 GeV (when 
the tree level expression is used). Our estimate of higher order effects, with the inclusion of the 
next-order effects in a s for massless quarks, shifts the curve for m b = 5 GeV in the direction of 
the 3 GeV result and amounts to about of 20% to 40% of the difference between the tree-level 
calculations with the two different masses. For both E and EM schemes we used the higher 
order results for the E scheme. 

For the JADE algorithm we have also plotted in fig. 3 the experimental results for R^ d 
obtained by the DELPHI group [2] on the basis of the data collected in 1990-1991. The 
experimental errors, due to the limited statistics analyzed, are rather large. However, one can 
already see the effect of the quark mass. If the 6-quark mass would be zero, one should obtain 
a ratio R 3 d constant and equal to 1. It is clearly seen from the figure that for y c < 0.08 the data 
are significantly below 1. For larger values of y c , the number of events decreases, the errors 
become too large and the data are consistent with 1 . When larger amount of data is analyzed 
and the experimental error is decreased, it will be very interesting to see if data will exhibit the 
different signs of the mass effect in R 3 d (positive for the E scheme and negative for the other 
schemes) as predicted by our parton level calculations (see fig. 3). 

In spite of the fact that the effect of the quark mass in R^ d has been seen, it is too early, 
in our opinion, to extract now the value of the 6-quark mass from the data. As discussed 
above the higher order corrections to R b 3 d are presumably rather large and should be included 
in the theoretical calculations. However, it is clear, that once the essential next-to-leading order 

5 With our choice of the normalization A^ 1 \y c ) = B{y c )/A, where B(y c ) is defined in [20]. 
6 For the EM algorithm this function has not yet been computed. To make an estimate of higher order corrections 
we will use in this case the results for the E algorithm. 
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Algorithm 


U&) 

y y y 


J2) 

'm 'm 'm 


EM 


-2.72 -14.64 -28.58 


-2.61 -13.54 -30.67 


JADE 


-2.01 -5.19 -13.25 


-1.90 -4.13 -15.42 


E 


4.68 19.04 25.97 


4.71 19.81 23.39 


DURHAM 


-1.69 - 4.76 -12.70 


-1.65 -4.28 -15.48 



Table 2: Results of the tree parameter fits of the functions By\(y c , r b )/A^°\y c ) = 
En=o k VA log" Vc in the range 0.01 < y c < 0.2 

corrections will be available and all LEP data will be included in the analysis, the ratios R^ d will 
certainly allow for a reasonable determination of the 6-quark mass and for a check of its running 
from m b to mz- 

To simplify the use of our results we present simple fits to the ratios By^ A (y c , r b )/A^°\y c ), 
which define R b 3 d at lowest order, for the different clustering algorithms. We use the following 
parameterization: 

B^ A (y c ,r b )/A(°\y c ) = £ k^log n y c , (3.11) 

n=0 

and the results of the fits for the range 0.01 < y c < 0.2 are presented in table 2. 

In fig. 4 we plot the ratios By^ A (y c , r b )/ A^°\y c ) as a function of y c for the different algorithms 
(dashed lines for m b = 5 GeV, dotted lines for m b = 3 GeV and solid curves for the result of 
our fits). As we see from the figure the remnant mass dependence in these ratios (in the range 
of masses we are interested in and in the range of y c we have considered) is rather small and 
for actual fits we used the average of the ratios for the two different masses. We see from these 
figures that such a simple three-parameter fit works reasonably well for all the algorithms. 

Concluding this section we would like to make the following remark. In this paper we 
discuss the zT-boson decay. In LEP experiments one studies the process e + e~ — >■ (Z*y*) — >■ bb 
and, apart from the resonant Z-exchange cross section, there are contributions from the pure 
7-exchange and from the 7 — ^-interference. The non-resonant 7-exchange contribution at the 
peak is less than 1% for muon production and in the case of 6-quark production there is an 
additional suppression factor Q\ = 1/9. In the vicinity of the Z-peak the interference is also 
suppressed because it is proportional to Q b (s — m|) {^/s is the e + e~ centre of mass energy). We 
will neglect these terms as they give negligible contributions compared with the uncertainties in 
higher order QCD corrections to the quantities we are considering. 

Obviously, QED initial-state radiation should be taken into account in the real analysis; the 
cross section for 6-pair production at the Z resonance can be written as 

= / <rMs')F(s'/s)ds' (3.12) 

where F(s'/s) is the well-known QED radiator for the total cross section [27] and, the Born 
cross section, neglecting pure 7 exchange contribution and the 7 — ^-interference, has the form 
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with obvious notation. Note that T b in this expression can be an inclusive width as well as some 
more exclusive quantity, which takes into account some kinematical restrictions on the final 
state. 



4 Angular distribution 

If the gluon jet can be identified with enough efficiency, an interesting quantity which is very 
sensitive to the IR behaviour of the amplitudes is the angular distribution with respect to the 
angle formed between one of the quark jets and the gluon jet 7 . If t?i (i? 2 ) is the angle between 
the quark (antiquark) and the gluon jets we define i? = min(i?i,i? 2 )- We want to obtain the 
angular distribution with respect to The starting point is eq. (3.2) where we change variables 
from one of the y 1 or the y 2 variables to To do this we take into account that the amplitude is 
completely symmetric in y 1 and y 2 , therefore we can restrict the integration only to the region 
y 2 > yi and add a factor 2. In that region i? = tV Therefore to obtain the distribution with 
respect to i? it is enough to obtain the distribution with respect to i?i but constraining the phase 
space integration to y 2 > yi . 

For y 2 > yi we can easily express yi in terms of cos i? = cos t?i as follows 

Vi (l-y 2 - cos^J(l - y 2 ) 2 - Ar b ) 

Vi = — J =-^ ■ (4.1) 

1 + y 2 + cosdJ(l - y 2 ) 2 - Ar b 



Adding the Jacobian of the transformation we find from eq. (3.2) (taking e = as this quantity 
is IR convergent) 



<*Tt r 2y 2 J(l-y 2 ) 2 -Ar b ) 

—f- = C b C g C F 2 / dy 2 e P s6 c e{y 2 - yi )sm$- ^ T A bg , (4.2) 

M J (1 + V2 + cos ^(1 - y 2 y - Ar b ) 

where y 1 is expressed in terms of cos ■& and y 2 using eq. (4.1). 

In order to see how large mass effects are in this angular distribution we define the following 
ratio of angular distributions: 

, , 1 dY\- I 1 dTt 

In the case of massless quarks the integration limits in eq. (4.2) can be found analytically 
for the JADE-type schemes and the result of the integration over y 2 is expressed in terms of 
logarithms involving i? and y c . We find 

1 dV% a 



vdM iMv*)> (4.4) 

where the function f^(y c ) is given analytically in appendix D for the JADE-type schemes and 
represented in fig. 5 for the JADE-type and the Durham algorithms for different values of y c 
(y c = 0.02 (solid line), y c = 0.04 (dashed line) y c = 0.06 (dotted line ) and y c = 0.08 (dash- 
dotted line)). We observe a very sharp peak, for both algorithms, in the region of 90°-100° 
depending on the value of y c , for y c = 0.04 the peak is at about 92° for the JADE-type algorithms 
and at about 99° for the Durham algorithm. We see that the absolute size of the peak is a factor 
two larger in the case of the JADE-type algorithms (for the same value of y c ) than in the case of 
the DURHAM scheme. This is due to the difference of phase spaces for two schemes. 



7 We thank J. Fuster for suggesting us the study of this observable. 
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For massive quarks, although the integrations can still be performed analytically in the EM 
scheme, some of the integration limits are solutions of polynomial equations of the third degree 
and the analytical result is not especially enlighting. Then, we have computed the ratio R b / by 
doing the one-dimensional integration in (4.2) numerically. 

Numerical results for R\f are presented in fig. 6 for the different algorithms for y c = 0.04 
and for both m b = 5 GeV (solid line) and for m b = 3 GeV (dashed line). In all cases we plot the 
ratios for the interval of angles for which the differential cross section is still sizable (see fig. 5), 
i.e. ■& 45° - 120° for JADE-type schemes and ■& 50° - 130° for the DURHAM algorithm. 
For small angles and m b = 5 GeV the effect can be as large as 10% of the ratio. Note, however, 
that the angular distribution, fig. 5, drops down rapidly for such small angles. In addition, since 
the ratio changes very fast in this region the exact size of the effect will depend on the angular 
resolution achieved at experiment. 

As in the case of ratios of three-jet event rates, R b 3 d , the variation of the ratio of angular 
distributions, R\f, for m b = 5 GeV and m b = 3 GeV gives a measure of the size of higher order 
corrections. 

We observe in all the ratios the irregular behaviour in the region where the massless angular 
distribution peaks. This is due to the fact that in the massive case the position of the peak is 
slightly shifted with respect to the massless case. The mismatch between the two peaks appears 
as a discontinuity in the ratio when seen from large scales. 

It will be interesting to see if data really follow these patterns for R\f. A preliminary analysis 
performedby the DELPHI group [28] seems to indicate that, indeed, data do follow these angular 
distributions, at least qualitatively, and exhibit the variations present in the different algorithms. 



5 Discussion and conclusions 

In this paper we have presented a theoretical study of quark-mass effects in the decay of the 
Z -bo son into bottom quarks. 

First, we have reproduced, with the complete mass dependences, the results for the inclusive 
decay rate of the Z — >■ bb + ■ ■ • to order a s by adding gluon bremsstrahlung from the 6-quarks 
to the one-loop corrected decay width of Z — >■ bb. Although the sum of the two contributions is 
finite, each of them is separately IR divergent. We used dimensional continuation to regularize 
the IR divergences and gave a complete analytical result in arbitrary space-time dimensions for 
each of the two contributions. 

The main contribution of this paper is, however, the analysis of some three-jet observables 
which are more sensitive to the value of the quark masses. 

For a slight modification of the JADE algorithm (the EM algorithm) we have calculated 
analytically the three-jet decay width of the Z-boson into 6-quarks as a function of the jet 
resolution parameter, y c , and the 6-quark mass. The answer is rather involved, but can be 
expressed in terms of elementary functions. Apart from the fact that these analytical calculations 
are interesting by themselves, they can also be used to test Monte Carlo simulations. For the 
EM, JADE, E and DURHAM clustering algorithms we have obtained the three-jet decay width 
by a simple two-dimensional numerical integration. Numerical and analytical results have been 
compared in the case of the EM scheme. 

We discussed quark-mass effects by considering the quantity 

pM it,fa.)/r 1 ' . m2 r , , 

R * =vm^ = + ^/ ( b ' Vc) 

which has many advantages from both the theoretical and the experimental point of views. In 
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particular, at lowest order, the function F(m b , y c ) is almost independent on the quark mass (for 
the small values of the mass in which we are interested in) and has absolute values ranging from 
10 to 35 (depending on y c and on the algorithm), where the larger values are obtained for y c of 
about 0.01. 

At the lowest order in a s we do not know what is the exact value of the quark mass that should 
be used in the above equation since the difference between the different definitions of the 6-quark 
mass, the pole mass, m b 5 GeV, or the running mass at the ra^-scale, fh b (m z ) ~ 3 GeV, 
is order a s . Therefore, we have presented all results for these two values of the mass and 
have interpreted the difference as an estimate of higher order corrections. Conversely one can 
keep the mass fixed and include in F(m b , y c ) higher order corrections already known for the 
massless case. According to these estimates the 0(a s ) corrections can be about 40% of the 
tree-level mass effect (depending on the clustering scheme), although we cannot exclude even 
larger corrections. 

By using the lowest order result we find that for moderate values of the resolution parameter, 
y c ~ 0.05, the mass effect in the ratio R^ d is about 4% if the pole mass value of the 6-quark, 
m b 5 GeV, is used, and the effect decreases to 2% if m b = 3 GeV. 

We have compared our predictions for R b 3 d for the JADE algorithm, with the results obtained 
from the 1990-1991 data by the DELPHI group [2]. Although the errors obtained in the analysis 
of this limited sample of data are rather large, especially for y c > 0.08, one clearly sees that 
for small values of y-cut (y c < 0.08) the experimental points are systematically below 1, thus 
clearly exhibiting the effect of the mass of the quark, as for massless quark R 3 d = 1. The 
size of the effect agrees roughly with the predictions. One can expect the reduction of the 
experimental error by, at least, a factor two when the data collected in 1992 are included in the 
analysis. Then, mass effects will be more clearly seen and it will be very interesting to see if 
data follow the different qualitative behaviour of the ratio R 3 d as a function of y c as predicted 
by the parton model calculations (positive effect for the E scheme and negative mass effect for 
the other algorithms). However, in order to extract a meaningful value of the 6-quark mass from 
the data it will be necessary to include next-to-leading order corrections since the leading mass 
effect we have calculated does not distinguish among the different definitions of the quark mass 
(pole mass, running mass at the m b scale or running mass at the mz scale). We believe that the 
future analysis of the whole LEP statistics and its comparison with the theoretical predictions 
for the three-jet ratios, which meet the future experimental precision, will allow for a good 
determination of the 6-quark mass at the highest energy scale and for a check of its running from 
m b to mz- 

The high precision achieved at LEP allows for a good separation of the gluonic and quark 
jets and a measurement of the angular distribution of the radiated gluon with respect to the quark 
momenta. This angular distribution has been calculated for massless quarks analytically for the 
JADE-type schemes and numerically for the DURHAM algorithm. We have studied the mass 
effects, for the different jet-clustering algorithms, in the quantity 



We have shown that, for a reasonable value of the resolution parameter, y c = 0.04, the mass 
effects in this ratio can be as large as 10% of the ratio for m b = 5 GeV (depending on the 
algorithm, the angle $ and the angular resolution). The larger values are obtained for small 
angles where, however, the angular distribution falls down very rapidly. A fit to this ratio can 
be used to extract the value of the 6-quark mass. It will be interesting to see if data really 
follow the predictions for the angular distributions and if the mass effects in the ratio of angular 
distributions are well described by our results. 




16 



Concluding, we have raised the question of the possibility of measuring the 6-quark mass 
at LEP by using three-jet observables. In our opinion, this is a big challenge for both experi- 
mentalists and theorists. Clearly, more work has to be done in order the precision of theoretical 
predictions meet the experimental accuracy, in particular order a 2 s calculations and studies of 
hadronization corrections including mass effects will be needed. However, this effort is worth 
since it will allow for an independent measurement of m b at much larger energies where, 
presumably, non-perturbative effects are negligible. 
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APPENDICES 



A Phase space in D = 4 — e dimensions 

The phase space for re-particles in the final state in D-dimensions [18] (D = 4 — 2e) has the 
following general form 



= l,n 



= (2*) D II r27r S 2E . ^ 2 -^ 2 ) Q (^)^ (?- E ^ • (a.2) 

i=l, n V / 1 \ i=l, n y 

Then doing several trivial integrations we have the following phase-space factor for the process 
Z ■ bb 

^ 2 " 4tt 2 r(2 - 26) { 4tt J ' (A3) 

where (3 = a/1 — 4r & with r & = va\jva 2 z , 

For the case of the decay into three particles, Z — >■ 66^, we have 

2 1 / 2 \ _2e 

d(PS 3 ) = mz XQ - ^ 0(h p )h- e d yi dy 2 , (A.4) 

V ; 16(2tt) 3 T(2 - 2e) V 4tt / V v) p y y ' V 7 

where the function Zip which gives a phase- space boundary in terms of variables 2/1 = 2(p 1 k)/ m 2 z 
and y 2 = 2(p 2 k)/m 2 z has the form 

K = 2/12/2(1 - 2/1 - 2/2) - ^(2/1 + 2/2) 2 • (A.5) 



B Inclusive decay rate functions 

In this section we collect the functions needed in section 2. The relevant form factors are: 



fAf 



C F |l + ^4Aog(c) 



2(3 



C 



+ 



2(3 

f e \og(r b ) + C F 
l + (3 2 ' 



log(c) , 
-2 



(B.l) 
(B.2) 



2 + (3 2 
2(3 



log(c) 



. Li 2 (c) + y " ^ log 2 (c) + log(c) log (1 - c) 



(B.3) 



In the expression for f^f, the first term, proportional to log(r&), comes because our election for 
the term proportional to the divergence. The vector form factor, f V f can be written in terms of 
the other two form factors, 

fv f = f Af + hf ■ (B.4) 
In terms of these form factors the functions F v and F A that appear in eq. (2.13) are 



Fr = p (<i^!) /v/ + §*,) , 

Fa = (3*f Af . 



(B.5) 
(B.6) 
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The functions that come from real bremsstrahlung can be written as follows, 

G v = (3 Q(3 - (3 2 )G P + G Vh ) , 



where 



G A = (3 (p 2 G P + G Ah ) , 



G P = Gp h + 2f £ (l + log((3)) 



(B.7) 
(B.8) 

(B.9) 



The terms proportional to f e come again from our choice of the coefficient of the divergence, 
and Gph is the finite part coming from the integration of the term proportional to h p in the 
amplitude 

Gp h = ^ [' C dzg(z)^ log (g(zfz(l + z)) . (B.10) 



2(3 

The result of the integration gives 



G 



ph 



C r 



-2 log 



4/3 3 \ 2 + (3\ 

+ 2 — log(c) 



1-/3 



l + (3 2 (I 



g ^log 2 (c) + y - Li 2 (c) - Li 2 (c 2 ) - 31og(c)log (1 + c) 



(B.ll) 



The functions Gy/, and G A u come from the integration of the h v and h A terms respectively 

(B.12) 



Cf f 1 ^ , ,. 2 ( 1\ C F / ^ 9-2/3 2 +/3 4 , , 



2/3 



and 



where 



G Ah = -(3 - /3 2 )GW + (1 - /3 2 )G^ , 
a O f 1/c , , x2 Of, , 2 3-2/3 2 -/3 4 1 N 

= ^ /< = Tv ^ — ^ — log(c) , 



(B.13) 
(B.14) 



C Three-jet event rate functions 

The functions Hy^ and \ which give the leading contribution to the three-jet decay rate in 
the EM algorithm, can be written in the following form 

H ( y\y c ,r b ) = C F J2 8i 
;=i,3 

;=i,3 

with 

9i = 9{Vc ~ Vc) 

2 = -e(y c -2r b )6(y c -y c ) 

#3 = 0{y c - y c ) 



(3-/3 2 ) 



Ks + Ky 



(3-/3 2 ) 



K v + (1 - (3 2 )K { A 



(C.l) 
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and y c — ^/rb)■ Here K l s corresponds to the soft part and K v ^ to the vector (axial) 

hard part. These functions are given by 



l-(3 2 



= Ay c (z^ - 1) - 2 \—^- + 2y c j log(z^) + 4/3 log ^_ ^ ? 



(C.2) 



+ 2(1 +/3 2 ) 



^ - U 2 (-z^) - Li 2 + U 2 (z^c) 



1 + 



y c 2 (i - ^ 2 ) 



32 



2(1 + zp)- 1 + 
1 + 



(l-/3 2 )(3 + /3 2 ) 



8 



8 ^ 

y 2 iog(^) 



(C.3) 



+ 



- + y 2 (l-^ 



(1-/3 2 ) 2 



16 



W - z ?) + C 1 + 



(l-/3 2 )(3 + /3 2 ) 
8 



log(^) 



(C4) 



if 



-2(l+/3 2 -2y c )log(z a ) + 4/31og 



1 - z a c 

2 



2 + log 



(l + /3 2 -2y c )z a (l + c) 2 



+ 2(1 +/3 2 ) 



log (z a ) log 



2(z a - c)(l - cz a ) 
(l+/3 2 -2y c )(l + c) 2 > 



Li 2 ( — ) + Li 2 (cz a ) 



(C.5) 



2\2 



Tf 2 



+ 



-(l + (3 2 -y c )y t 

(l-/3 2 )(3 + /3 2 ) 
8 

/ (l-/3 2 )(7 + /3 2 ) 



(1 ~ z a ) , (I-/? 2 ) 



+ 



32 



■0 



a a.) 



(2-z- 1 + z a -4(l + z a )- 1 ) 



+ yc + (3 2 y c -y 2 log(z a ) 



(C.6) 



K 2 



+ 



(1-/3 2 ) 2 _x . (l+/3 2 -y c )y c (l-z a ) 
[z a - z a ) + 



16 

(1 - (3 2 )(3 + (3 2 ) 
8 



2 (1 + z a ) 

(-l + 2(l + z a )- 1 +log(z ce )) 



(C.7) 
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(1+gV (1 - /3 2 - 4/3 2 z 7 ) (1 - z 7 ) (1 + 3/3 2 - 2(1 - /3 2 )z 7 ) 



+ (1+/3 2 ) 



log 2 



z 7 (1 + 2z 7 ) 

.1 + Z 7 ) +2Ll2 (i + ^ 7 , 



(1 + 2z 7 ) 



(C.8) 



8 



(1 + 2z 7 ) 2 (1 + 2z 7 ) 



log(z 7 ) 



(C.9) 



(1+/3 2 ) 2 / l-z 7 V 
^ " 8 l.l + 2zj ' 

where we used the following notation, 



^- (l - y c - 2r b - y/(l - y c f - 4r 6 ) 



(CIO) 



l-2r b - 2y c ' 
1-/3 



In the limit of massless quarks, r b — >■ 0, from the functions fT-^ ^ and given above we 
obtain 

H$\y e ,r b - 0) = #i 0) (y c ,7- 6 - 0) - A<°> , 
Here the function A^°\y c ) is the known result [20,21] for the JADE algorithm 



r(°) 



(Cll) 



(C.12) 



A^\y c ) 



2C t 



+ 21oe 



y+2-6 yc -- yc + (3-6 yc )log Vi _ 2yc 



+ 4Li 2 



i - Vc 



(C.13) 



The function A(y c ) given in refs. [20,21] differs from our A^°\y c ) in a factor 2 because we 
chose a different normalization for it. 



D Angular distribution functions 

The angular distribution studied in section 4 is given, in the massless case, by the function 
f#(y c )- In the JADE-type algorithms it can be written as follows 



fo(y c ) = C FS m(#) £ eifi(y c ) 
;=i,2 



(D.l) 



where the theta; functions have the form 



0l = \co S ^ + 77 ^\ e( {1 ~ 6y < \P - cos A , 



(i-yc). 



(i + z/c) 2 



61, = 61 



(i-y c ) 



cos ■& 6 1 cos i? + 



1 ~ Vc 

l+Vc 



21 



and 



h{y c ) 



In these equations we defined 

6 ■- 

X X : 

x 2 ■- 



(1 + &) 2 



y c + 86 - 3 



log Q +(l + 6 + 26 2 )log 



(1 - y c ) 2 - Aby, 

' b + xA 
b + x 2 



(1 + &) 2 



{Vc - b){y 2 c + 2y c - 2by c - 5b 2 ) 



6(1 + 2b) log 



1 + COS 1? 



2b 



b + y c 



4(6 + y c ) 
- log 



1-6 

.1 - y c 



^ (i - y c - V(i-y c ) 2 -46y c ) , 

^ (l - y c + ^(1 - y c ) 2 - 46y c ) . 



(D.2) 



(D.3) 
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Figure 2: The phase space for Z — >■ bbg in the plane yi and y 2 with cuts (y c = 0.04 and 
y c = 0.14) for the different algorithms. The mass of the quark has been set to 10 GeV to 
enhance mass effects in the plot. 



26 




Figure 3: The ratios (see eq. (3.8)) for the four algorithms. Solid lines correspond to 
m b = 5 GeV, dashed lines correspond to m b = 3 GeV and dotted lines give our estimate of 
higher order corrections to the m b = 5 GeV curve. For the JADE algorithm we have also 
included the results of the analysis of the data collected during 1990-1991 by the DELPHI 
group [2]. 
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Figure 4: The functions By /A^ and B^ /A^ for the four algorithms. Dashed lines for 
m b = 3 GeV, dotted lines for m b = 5 GeV and solid lines for our three-parameter fit, eq. (3.1 1). 
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Angular distribution (m =0) 
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Figure 5: Normalized angular distributions (eq. (4.4)) with respect to the angle formed between 
the quark and the gluon jets for the massless case for JADE-type and DURHAM algorithms. 
Solid line for y c = 0.02, dashed line for y c = 0.04, dotted line for y c = 0.06 and dash-dotted 
line for y c = 0.08 
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Figure 6: The ratios of angular distributions (see eq. (4.3)) for y c = 0.04 for the different 
algorithms. Solid line for m b = 5 GeV and dashed line for m b = 3 GeV. 
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